Quantum Transport 

Lecture Notes 

Yuri M. Galperin 

Lund University 
November 1998 



E-mail: yurig@teorfys.lu.se 
Permanent address: 



Department of Physics, P.O. Box 1048 Blindern, 0316 Oslo 
Phone: +47 22 85 64 95, E-mail: iouri.galperine@fys.uio.no 



ii 



Preface 



This course is an attempt to provide an overview of basic concepts of quantum transport 
through modern small-size structures. It is a very hot topic at present time because it is 
relevant to fundamental principles of quantum mechanics and statistical physics, as well 
as to various applications in modern electronics. 

Several excellent books [II, i, 0, 0] and review articles (e. g. [0, H, 0, QU, O, 03, 
□3, [Q], H3, HB, IT3, HSI, H3]) are written on this subject. This list is far from being complete. 

In the present course I try to select only few topics from the broad area to provide a 
general introduction to the subject. The reference list includes only few selected papers 
rather than a more or less complete overview of the literature. I tried to concentrate on 
the basic concepts rather than on historical aspects. 

Any comments and criticism will be gratefully appreciated. 



October 1998 
Yuri Galperin 
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Chapter 1 

Preliminary Concepts 



1.1 Two- Dimensional Electron Gas 

An important system where quantum effects were observed is two-dimensional electron gas 
(2DEG). There are two basic systems where 2DEG has been studied. One of them is Si 
MOSFETs (metal-oxide-semiconductor field-effect transistors). A very good review of such 
systems is given in Ref. [□]. A typical device is shown in Fig. [T7T| . A (100)Si surface serves 
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Figure 1.1: Band diagram showing conductance band Ec, valence band Ey and quasi- 
Fermi level Ep. A 2DEG is formed at the interface between the oxide (SiO^) and p-type 
silicon substrate as a consequence of the gate voltage V g . 

as a substrate while SiC>2 layer behaves as an insulator. 2DEG is induced electrostatically 
by application a positive voltage V g . The sheet density of 2DEG can be described as 



71, 



<Vg-V t ) 



where V t is the threshold voltage for the barrier's creation 

Another important systems with 2DEG involve modulation-doped GaAs-AlGaAs het- 
erostructures. The bandgap in AlGaAs is wider than in GaAs. By variation of doping 
it is possible to move the Fermi level inside the forbidden gap. When the materials are 
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put together, a unified level of chemical potential is established, and an inversion layer is 
formed at the interface. 
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Figure 1.2: Band structure of the interface between n-AlGa As and intrinsic GaAs, (a) 
before and (b) after the charge transfer. 



The 2DEG created by a modulation doping can be squeezed into narrow channels 
by selective depletion in spatially separated regions. The simplest lateral confinement 
technique is to create split metallic gates in a way shown in Fig. A typical nanostructure 
is shown in Fig. |1.4| . 



1.2 Basic Properties of Low-Dimensional Systems 
Wave Functions 

Let us direct 2-axis perpendicular to the plane of 2DEG. The wave function can be decou- 
pled as 

where r is the vector in plane of 2DEG. Throughout our considerations we will assume that 
all the distances are much larger than interatomic distance and thus we will use the effective 
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Figure 1.3: On the formation of a narrow channel by a split gate. 



Figure 1.4: Scanning electron microphotographs of nanostructures in GaAs-AlGaAs het- 
erostructures. Taken from M. L. Roukes et ai, Phys. Rev. Lett. 59, 3011 (1987). 

mass approximation. A good approximation for the confining potential is a triangular one, 

oo at z < 



U(z) 



i 

Fz at z > . 



Then one can write the Schrodinger equation for the wave function x( z ) as 

W + ^ E - FZ)X = ' (L1) 
Instead z we introduce a dimensionless variable 

E\ {2mF\ 1/3 



C=[z , 



\ (2mF\ 
) \J^J 



The quantity 

plays the role of characteristic localization length in z direction. Then Eq. (|1.1|) acquires 
the form 

x" - Cx = o 
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which should be solved with the boundary conditions of finiteness at infinity and zero at 
z = 0. Such a solution has the form 

x(C) = AAi(C). 

Here Ai(£) is the Airy function defined as 

1 f 00 

Ai(C) = —7= / cos(m 3 /3 + u() du . 
V 71 " Jo 

For large positive ( it decays exponentially, 

Ai(C) - 2^e- (2/8, <"' a , 
while for large negative zeta is is oscillatory, 



Ai «~K^ sil # |3/2 + £) 



The energy spectrum is defined by the roots ( n of the equation 

Ai(C) = , -> £ n = -EoCn • 

Here 

£ = 



#2^2 \ 1/3 



2m 

We have £i ~ —2.337, C2 ~ —4.088. The normalization constants t4„ for each level are 
defined as 

A' 1 = / dz\ X n(z)\ 2 . 



Normalized electron densities y4 n |x„(2;)| 2 are shown in Fig. 1.5 . Each level creates a sub- 
band for the in-plane motion, the energy being 

h 2 k 2 

E n k = E n + E(k) = E n + — — . 

2m 

Note that the effective mass m is considerably smaller than the mass of a free electron. 



Density of States 

The density of states g(e) is defined as number of states per the energy interval e, e + de. 
It is clear that 

a 

where a is the set of quantum numbers characterizing the states. In the present case it 
includes the subband quantum number n, spin quantum number a, valley quantum number 
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v (for n-type materials), and in-plane quasimomentum k. If the spectrum is degenerate 
with respect to spin and valleys one can define the spin degeneracy v s and valley degeneracy 
v v to get 



n J 



(27T 



Here we calculate the number on states per unit volume, d being the dimension of the 
space. For 2D case we obtain easily 

u a u v m ^ 



■i i 



Within a given subband it appears energy-independent. Since there can exist several 
subbands in the confining potential (see Fig. |T76| , inset), the total density of states can be 



represented as a set of steps, as shown in Fig. |L6| . At low temperature (kT <C Ep) all the 
states are filled up to the Fermi level. Because of energy-independent density of states the 
sheet electron density is linear in the Fermi energy, 

f v s v v mE F 

n s = M - — h const 

znn 2 
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^1 ^2 Ep E3 Energy 

Figure 1.6: Density of states for a quasi-2D system, 
while the Fermi momentum in each subband can be determined as 



k Fn = i ^2m(E F - E n ) . 

Here M is the number of transverse modes having the edges E n below the Fermi energy. 
The situation is more complicated if the gas is confined into a narrow channel, say, along 
2/-axis. In a similar way, the in-plane wave function can be decoupled as a product 

where N is a proper normalization factor, the energy being 



E 



n,s,k 



E n + E s {k x ) =E n + E s + 



h 2 k 2 x 

2m 



Here E ns = E n + E s characterizes the energy level in the potential confined in both (z and 
y) directions. For square-box confinement the terms are 

= (snhf 
s 2mW 2 ' 

where W is the channel width, while for the parabolic confinement U(y) = (l/2)muj 2 y 2 
(typical for split-gate structures) 

E s = (s- l/2)hoj . 

It is conventional to introduce partial densities of states for the states with k x > and 
k x < 0, g ± , respectively. We have, 



07 W 



u 8 v v ( dE s {k x 



2tt 



v s v vy /m 



The total density of states is 



2 VHh V^EZ ' 



E 



;i.2) 



;i.3) 



The energy dependence of the density of states for the case of parabolic confinement is 
shown in Fig. |1.7| . 
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Figure 1.7: Density of states for a quasi-lD system (solid line) and the number of states 
(dashed lines). 

Motion in a perpendicular magnetic field 

2DEG in a perpendicular magnetic field gives an example of O-dimensional electronic sys- 
tem. Indeed, according to the classical theory the Hamilton's function of a charged particle 
in an external electromagnetic field is 



where is the scalar and A is the vector potential of the field, and p is the generalized 
momentum of the particle. According to the rules of quantum mechanics, one should 
replace the canonical momentum p by the operator 



and add also an extra spin term — fiH where fi = /igs/s. Here fis = e/2mc is the Bohr 
magneton while s is the spin operator. Generally, interaction with periodic potential of the 
crystalline lattice leads to renormalization of the spin splitting \i B — > fi = gffiB where gf is 
called the spectroscopic spin splitting factor. 
Finally we get, 




p — > p = —iftV 



n 





— (A. P + P .A) + 



I - ^s-H + e0. 

- s 
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Since 

p ■ A — A ■ p = —ih divA , 
those operator commute if divA = 0. It holds in a uniform field with 

A = -H x r. 

2 

The wave function in a magnetic field is not uniquely defined: it is defined only within the 

gauge transform 

A^A + V/, 0-^0--^, 

c at 

where / is an arbitrary function of coordinates and time. Under such a transform only the 
phase of wave function is changed by the quantity ef /he that does not affect the observable 
quantities. 

In classical mechanics, the generalized momentum of the particle is related to its velocity 
by the Hamilton equations, 

mv = p — eA/c . 

According to the quantum mechanics we arrive at a similar expression. However different 
components of velocity do not commute, the commutation rules being 

{v x ,v y } = i(eh/m 2 c)H z , 
{v y ,v z } = i(eh/m 2 c)H x , 
{v z ,v x } = i(eh/m 2 c)H y . 

That means that the particle cannot simultaneously have definite velocities in all three 
directions. 

Let us determine the energy levels in a 3-dimensional system embedded into a uniform 
magnetic film with a vector potential 

A x = -Hy , A y = A z = . 

The Hamiltonian then becomes 

n = J_ U + e l£\ + El + &_ _ a S H . 

2m \ x c J 2m 2m s 2 

First, the operator s z commutes with the Hamiltonian. Thus ^-component of spin is 
conserved and can be replaced by its eigenvalue a. Thus one can analyze the Schrodinger 
equation for an ordinary coordinate function, 

/j, 

ip aHtp = Eijj . 

s 

It is naturally to search for solution in the form 

i) = e i{PxX+PzZ),h <p{y) . 



1 

2m 



('II \ ; 

Px + —y ) +p y + p z 
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The eigenvalues p x and p z take all values from — oo to oo. Since A z = we get 

p z = mv z . 

Thus the motion along magnetic field in 3D system is not quantized. For a motion in the 
xy-plane we have the following Schrodinger equation, 



2m 



E+^H 

s 



K_ 
2m 



~-muj 2 c (y - y )' 



Here 



Vo 



-cp x /eH = -a 2 H k x , a H = {ch/eH) 



1/2 



0. 



\e\H/mc . 



;i.4) 
;i.s) 



Since this equation is the same as the Schrodinger equation for a harmonic oscillator, we 
obtain 



E = (n + l/2)huj c - (imj/s)H + p 2 z /2m. 



n 



0,1,... 



;i.6) 



The first term gives discrete levels which corresponds to the finite motion in the xy-plane, 
they are called Landau levels. For an electron, fi/s = —\e\h/mc, and the energy spectrum 
reads as 



E 



n H h a huj c + 



2m 



The eigenfunctions (p n (y) are 

My) z 



exp 



7T 



(y-yo) 2 ' 




'y-yo' 


L 2a|r . 







;i.7) 



;i.8) 



Here H n is the Hermite polynomial. 

In classical mechanics the motion in a magnetic field in xy-plane takes place in a circle 
about a fixed center. Here the conserved quantity y corresponds to y coordinate of the 
center of the circle. It is easy to see that the combination 

xq = cpy/eH + x 

is also conserved, it commutes with the Hamiltonian. The quantity Xq corresponds to a clas- 
sical x coordinate of the circle center. However, the operators and xq do not commute. 
That means that the coordinates xq and yo cannot take definite values simultaneously. Q 
One can ask: why the coordinates x and y are not equivalent? The reason is that 
the wave functions ( |1.8| ) correspond to the energy independent of k y . Consequently, any 
function of the type 



1 In a classical mechanics, the radius of the circle is r c = cmv t /eH = v t /uj c . Here v t is the tangential 
component of the velocity. Thus we have, 



y = Vo + r c (v x /v t ) , x = x - r c (v y /v t ) . 
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corresponds to the same energy and one can chose convenient linear combinations to get 
correct asymptotic behavior. 

To calculate the density of states in a magnetic field first we should count the number 
of the values k y corresponding to the energy e a (the so-called degeneracy factor). As usual, 
we apply cyclic boundary conditions along y and z -axes and get 

2tt _ 2tt 

K x — Tly, K z — 7l z . 

At the same time, we assume that the solution exists only in the region 

< y < L y . 

So, the degeneracy factor is 

^ X \U |max _ L x x _ L y L x 

2^ ~27ral y ° ~ 2va? H 1 9j 

This is very important relation which shows that one can imagine Landau states as cells 
with the We will come back to this property later. 

Now it is easy to calculate number of states in a 3D system treating the k z variable as 
for the usual ID motion 

2\k z \L z 2\[2~mL z / — — — — — 

for each state with a given N. Finally, the total number of states per volume for a given 
spin is 

Z s (e) = Z s n{e) = jp^j- Ve - hco c (N + 1/2) 
at y Z7T ) aa H N 

where one has to sum over all the values of iV with non-negative e — t\oj c (N +1/2). The total 
number of sates is Z(e) = 2Z s (e). To get DOS one should to differentiate this equation 
with respect to e. The result is 

dZ(e) V2^n~ ^ 6[e - hw c (N + 1/2)] 
9s[£) (2n)W H ^ y/e-tUj c {N + l/W 

Here 

{1 for x > ; 
1/2 for x = 0; 
for i<0 

is the Heaviside step function. To take the spin into account one should add the spin 
splitting ±HB9fH to the energy levels. If we ignore the spin splitting we can assume 
spin degeneracy and multiply all the formulas by the factor 2. We take it into account 
automatically using g(e) = 2g s (e). 
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Figure 1.8: Landau levels as functions of p z (left panel) and of H (right panel). The Fermi 
level is assumed to be fixed by external conditions. 

The behavior of the density of states could be interpreted qualitatively in the following 
way. The Landau levels as functions of magnetic field for a given value of p z are shown in 
Fig. |1.8| . As a function of magnetic field, they form the so-called Landau fan. The Fermi 
level is also shown. At low magnetic fields its dependence on magnetic field is very weak. 
We see that if magnetic field is small many levels are filled. Let us start with some value 
of magnetic field and follow the upper filled level N. As the field increases, the slopes of 
the "fan" also increase and at a given threshold value for which 



As the field increases the electrons are transferred from the iV-th Landau level to the other 
ones. Then, for the field H N _i determined from the equation £jv-i(-ffjv-i) = € f the (N—l) 
becomes empty. We get 



Here m c is the so-called cyclotron effective mass which in the case of isotropic spectrum is 
the same as the density-of-states effective mass. We observe that DOS in a given magnetic 
field oscillated with the increase in energy just similar to the case of quasi ID systems. 
Here Landau sub-bands play the same role as the modes of transverse quantization for 
quantum channels. 

For a 2DEG the motion along z-direction is quantized, and instead of e %VzZ l n we have 
Xs{z). The means that for each subband of spatial quantization we have a sharp Landau 
level, the density of states (per area) being 



Thus the density of states has sharp maxima at the energy levels that is a feature of 
so-called O-dimensional system. In real samples the peaks are smeared by disorder. 



£n(Hn) — &f- 



eh N 



so 
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1.3 Degenerate and non-degenerate electron gas 

At equilibrium the states are filled according to the Fermi function 

fo(e) ' 



exp[(e - fi)/kT\ + 1 ' 

where fi is the chemical potential while k is the Boltzmann constant. The chemical potential 
is determined by the normalization to the total number of electrons as 

poo 

n = / g(e) /o(e) de 
Jo 

where n is the electron density. At zero temperature the chemical potential is called the 
Fermi energy, ep- The graph of the Fermi function and its energy derivative is given in 
Fig-HJ 




Figure 1.9: The Fermi distribution (solid line) and its energy derivative multiplied by kT 
(dashed line) for (/kT = 10. 

Since at T = 

/o(e) « 9(e - C) , 
the Fermi energy is given by the equation 

rep 

n= g{e)de. (1.10) 
Jo 

The limiting case T = is actually means the the inequality kT <C ep is met. In the 
opposite limiting case, kT ^> e^, we get 

POO 

/o(e) w e^ kT , n = e^ kT / g(e) e^ kT de . 

Jo 
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Thus, 

1 1 r°° 

/o(e) = A(T) e^ kT , = - J g(e) e^ kT de . (1.11) 

This distribution is called the Boltzmann one. 

1.4 Relevant length scales 

One can discriminate between several important length scales in low-dimensional systems. 
They are shown in the Table |1.1| . 



1 mm 


Mean free path in the quantum Hall regime 


100 jum 


Mean free path/Phase relaxation length 
in high-mobility semiconductor at T < 4 K 


10 //m 




1 fim 


Commercial semiconductor devices (1990) 


100 nm 


de Broglie wave length in semiconductors. 
Mean free path in polycrystalline metallic films 


10 nm 




1 nm 


de Broglie wave length in metals 
Distance between atoms 


1 A 





Table 1.1: A few relevant length scales. Note that 1 /im = 10 6 m = 10 4 cm; 1 nm = 
10~ 9 m = 10 A. 

The above mentioned scales have the following physical meaning: 
De Broglie wave length, A. This length is defined as 

p k 

where p (k) is the typical electron momentum (wave vector). For Fermi gas the 
characteristic momentum is just the Fermi momentum. For the case of a single filled 
band in 2DEG, 

A = 2n/kF = \/2nJn~ s 
where n s is the sheet density. For the Boltzmann gas, p ~ \f2mkT, and 
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Mean free path, I. This is a characteristic length between the collisions with impurities 
or phonons. It is defined as 

i = VT tI 

where v is the typical velocity while r tr is the so-called transport relaxation time. It 
is defined as 



— oc / d0sin0W(0)(l-cosi 

Ttr y 



where 6 is the scattering angle while W(0) is he scattering probability. Usually the 
transport is characterized by the mobility 



er tI 

u = 

m 



The physical meaning of mobility is that a typical electron drift velocity acquired in 
an external electric field E is given by the relation 

v d = uE . 



Phase-relaxation length, L v . This is a specially quantum mechanical relaxation length 
which has no analogs in classical physics. Namely, classical motion can be described 
as evolution of the probability to find a particle at a given point at a given time. 
However, in quantum mechanics the state is characterized by the wave function which 
has a phase. The phase is important in the so-called interference phenomena, where 
the electron wave functions having different pre-history are collected at the same 
point. If the phases of the waves are not destroyed, a specific quantum interference 
phenomena can be observed and important. The phase-relaxation time, t v , describes 
relaxation of the phase memory. 

It is clear that scattering against any static spin-independent potential cannot lead 
to the phase relaxation. Indeed, in any stationary potential the equations of motion 
are time-reversible. The only processes which can be responsible for phase relaxation 
are the ones which broke the symmetry with respect to time-reversal. Among them 
are inelastic scattering by phonons, electron-electron collisions, spin-flip processes, 
etc. An important feature of such processes is that an electron suffers many elastic 
collisions during a typical time r v . Since it moves diffusively a proper way to estimate 
the relevant length L v is as follows: 

Lip = a/ Dt v , 

where D = (l/d)v£ is the diffusion constant (d is the dimensionality of the electron 
gas). 

Thermal dephasing length, L T . The above mentioned relaxation process is relevant to 
the interference of the wave functions belonging to a single-electron state. However, 
interference can be also important for the interaction of two electrons having close 
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energies. Indeed, if the energy difference between the electrons is ~ kT they travel 
almost coherently during the time H/kT. Thus the characteristic length of coherent 

propagation is estimated as 

L T = ^JhD/kT. 

Comparing mean free path £ with characteristic dimensions of the system, L, one can 
discriminate between diffusive, f C L and ballistic, £ > L, transport. Such a classification 
appears incomplete in the situation where different dimensions of the sample are substan- 
tially different. The situation is illustrated in Fig. |1.1U| for the case where the length L of 
the sample is much larger than its width, W. If phase coherence is taken into account, the 




Figure 1.10: Electron trajectories for the diffusive (£ < W,L), quasi-ballistic (W < £ < L) 
and ballistic (£ > W, L) transport regimes. From [EH]. 

scales L v and L? become important, and the situation appears more rich and interesting. 
Mesoscopic conductors are usually fabricated by patterning a planar conductor that has one 
very small dimension to start with. Although some of the pioneering experiments in this 
field were performed using metallic conductors, most of the recent work has been based on 
the gallium arsenide (GaAs)-aluminum gallium arsenide (AlGaAs) material system. Some 
important parameters for such systems are shown in Fig. |1 . 1 1| . 
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Figure 1.11: Electronic properties of the 2DEG in GaAs-AlGaAs and Si inversion layers. 
From IDHI. 



Chapter 2 
Diffusive transport 



2.1 Classical transport in diffusive regime 
Boltzmann equation 

According to classical physics, one can specify the coordinate r and the momentum p of 
the particle. Thus one can calculate the non-equilibrium distribution function, / p (r, t), the 
current density being 



j(r,f) = e J(dp)vf p (r,t). 



Here we denote (dp) = d d p/(2nh) d . The distribution function is defined by the Boltzmann 
equation 

df p (r,t) _ df dfdr dfdp 

+ ~ 77 + 7^ 7. r -t C oll 



dt dt dr dt dp dt 

df df df T 

at or op 
= (2.1) 



Here 



F = e (e + i[v x H] 



is the force acting upon the electrons, v = de p /dp is the (group) electron velocity, while 
J coll is the collision operator. It expresses the changes in the state due to quantum collisions 
and can be expressed through the the transition probability Wif between the initial state 
(i) and the final state (/), 

Icon(fa) = i W ^fa(l - fa') ~ W^f^l ~ f a )} ■ (2.2) 

a' 

Here a denotes the electronic state (in our case a = {p, s}). 

17 



18 



CHAPTER 2. DIFFUSIVE TRANSPORT 



In a stationary situation and in the absence of external fields the solution of the Boltz- 
mann equation is the Fermi function. Substituting it into the collision operator we obtain 
a very important relation between the probabilities of the direct and reverse processes, 

W aa ,e- e «' lkT = W Q , a e- ea/kT . 



Drude formula 

The simplest problem is to calculate the linear response to a small stationary electric field, 
E. Since we assume E to be small, it is reasonable to search solution as 

/ p (r) = / (ep) + , |/ (1) |«/ . 
Since I coU (fo) = we come to the integral equation 



The linearized collision operator has the simplest form for elastic processes when W pp i = 
W p 'p. Then, 

w/ (1) )=E w w(tf ) -^ 1) ) • 
p' 

Since we are interested in the function odd in p (to create a current) and it must be a 
scalar it is natural to assume that 

/Woe (p- u), v = E/E . 

Under such an assumption 

T tr r tr , \ Pv J 

p 

The quantity r tr is called the transport relaxation time. If the material is isotropic then W 
is dependent only on the scattering angle 6 between p and p'. From the simple geometrical 
construction shown in Fig. |TT] we observe that 

Pu = p' cos ^ cos ^5 Vv = P cos (j> ■ 
Consequently, we can express the relaxation time in the form 

1 , A 

Using the expression for the relaxation time we can write 



9( e p)l [ d6 sm6 (l-cos6) W(6). (2.3) 
2 Jo 



/W = r tr e(E • v) (-^) - j = e 2 J (dp) v(E ■ v) r tr (e p ) (- 



dfo(e P Y 
de r 



p 



2.1. CLASSICAL TRANSPORT IN DIFFUSIVE REGIME 
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Figure 2.1: On the calculation of the transport relaxation time. 
As a result, we arrive at the Ohm's law, j = crE with 

* = [ W)(^ r > £ (-^) • (2.4) 

Here means the projection of the velocity on the direction of the eelctric field, while 
(. . .) e means the average over the surface of a constant energy e, 

(A(p)) e - ^ (rfp)A(p)<5(e_ep) 



f(dp)6(e-e p ) ■ 
The quantity 

D = «r tr ) e = ^v\ x = l -vi (2.5) 

has an explicit physical meaning. This is just the diffusion constant for the electrons with 
a given energy e. Indeed, for low temperatures we get 

a = e 2 g{t F )D{t F ) . (2.6) 

On the other hand, the phenomenological expression for the current density in the presence 
of the density gradient reads as 

j = <rE - eDVn , 

where D is th diffusion constant. In the equilibrium the current is zero and 

eDVn = eDg(e F )V( = -aVcp . 

At the same time, the electro-chemical potential £ + ip/e, must be constant and Vip = 
— eVC- Substituting Eq. (|2.6|) we identify the quantity D(e F ) with the diffusion constant. 
Eq. ( |2.6|) is known as the Einstein relation. 



20 



CHAPTER 2. DIFFUSIVE TRANSPORT 



2.2 Linear response in quantum mechanics 

Let us find a linear response to the perturbation Q 

= lim Ae-^Ae"* . (2.7) 

17— >0 

According to time-dependent perturbation theory, the change in any quantity (B) in the 
Heisenberg representation can be written as 

5(B(t)} = l -£jt'(g\[H%(t'),B H (t')}\g} 

f dt'Y,{{9\n^\n){n\B H \g) - (g\B H \n)(n\H^\g)} . (2.8) 

J-oo 



I 

h 



Here \n) means a complete set of eigenstates for the unperturbed Hamiltonian TC, g stands 
for the ground state, while subscript H stands for Heisenberg representation with respect 
to non-perturbed Hamiltonian, TC. Going to the Schrodinger representation as, 

(g\A H \n) = e^ nt (g\A\n) = e^ nt A gn , ftoo gn = E g - E n 

we find that 5B(t) = 5B u e~ iwt with 

^ = l i m^E(- <gH " >( " |B|9> + <9|B| " >< " |a|g> |- (2-8) 

V^O ft / — ' ^ UJ + UJ ng + 77] UJ — UJ ng + 17] ) 

Defining linear response as SB^ = Xba{u)X we find 

= Um n \ E I" t 9nB l- + B9nA l } ■ ^ 
V^O ft y UJ + UJ ng + 17] UJ — UJ ng + 17] ) 

Let us use this general formula for the response of a dipole moment density er/V in the 
direction j3 to the electric field E directed in the direction a. Here V is the volume of the 
sample. Since the perturbation Hamiltonian is eE • r we find 

X/3a = ^lim^N ^L_ + ^ . (2.11) 

ftV r?->0 I + Co> n5 + Zr/ Co" - CU n9 + 2?7 J 

Since the complex dielectric function is 

1 + Aixx(u) = 1 + 47r«T(cj)/u; 

we obtain 

= ioo^- lim V «" "« • (2.12) 

V/i ^0 ^ UJ + UJ na + 17] UJ — UJ na + 17] 



UJ + LU ng +17] UJ — LU ng + 17] 

1 Here I partly follow the derivation given in Ref. [01]. 



2.2. LINEAR RESPONSE IN QUANTUM MECHANICS 
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In particular, 

2 

= E W" + - - ^>fnCj (2-13) 

n 

Now we can generalize the presentation for finite temperature doing inverse Fourier trans- 
form. Since S(u) = J dt exp(— iuot) we can write 

->2 r°° 



Ua Pa {u) = ^J™Jt(g\[r a (0)y(t)]\gy 

dt([r a (0),r^(t)]) r e-^. (2.14) 



cue 2 rco 



2Vh 

Here we have replaced ground state average by the thermal one to alow for finite temper- 
atures. 



Fluctuation-dissipation theorem 

Our derivation allows to obtain an important fluctuation- dissipation theorem relating fluc- 
tuations in the system with the dissipation. Let us characterize fluctuations by the factor 

/oo 
dt(r«(0y(t)) T e^ = ^iZ&Sfa - Un9 ) P g 
00 g,n 

where P g is the thermal weight of the state g. Using equilibrium thermal weights^ for Fermi 
particles, P g = fo(E g )[l — f (E n )}, where fo(E) is the equilibrium Fermi distribution, we 
obtain 

S^H = e ^Sp a {u) [I - e-^ kT ] (2.15) 

This theorem is often used in some other form which can be obtained using the relations 
between the operators v = r. Assuming hu <C kT we get 

2 roo 

^ axx{uj) = 2VkfJ rf ^(0K(t)>Te^ (2-16) 
Writing (v x (0)v x (t)) T as a Fourier transform of velocity fluctuations, S v (u), we get 

e 2 S v {u) kT 



V 7T 



which is just the famous Nyquist- Johnson relationship between the current noise spectrum 
and the conductivity. 

2 Here we actually assume the transport through the device to be incoherent. 
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At finite temperature using the expression for P g we get 



e S v (u) Hlu ( 1 1 



Ka(cj) (2.17) 



V 7T V 2 - 1 

Finally, since the contribution of an electron to the current is ev/L, 

Si(u>) = — coth (%) 

v ; 2tt V 2 ^y v ; 

Here G(cu) = aA/L is the conductance of the system (.A is the area, while L is its length). 
Derivation of Drude formula 

Now we can derive the Drude formula explicitly. Coming back to the expression fl2.14J) and 
introducing velocity operators we get 

&<Tap(w) = y / dEg(E) h{E) - f ^ E + M dt e -*-(^(0y(t)) B . (2.18) 

Here g(E) is the density of states while subscript E means that the average is calculated 
over all the states with given energy E. 
The quantity 

D aP {E) = - / dte-^(v a (0y(t)) E (2.19) 

is nothing else that diffusion constant. Indeed, assume that we have a small gradient in 
electron density, n(x) = n + ex. Then the particle current is 

j x = Hm(^(t = 0)n[x(t = -At)]) = Hm c{v x (0)x(-At)) 

/•At 

= -c Hm / dt (v x (0)v x (-t)} . 
A '^° Jo 

In this way we obtain the classical expression for diffusion constant. Finally, at hu <C \l 
we arrive at the formula 



9f^/3 = e 2 J dEg(E) D a0 {E) . 



Substituting 

(v 2 ) r°° 

D aP {E) = ^ / e -^*- |t|/Tt 

2 J-oo 



V a) f°° „— tort— |t|/Ttr _ ^ T tr 

d 1 + (cur tr ) 2 



(where ti is dimensionality of the electron motion) we finally arrive at the Drude formula. 



2.3. SHUBNIKOV-DE HAAS EFFECT 
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2.3 Linear response in magnetic field. Shubnikov-de 
Haas effect 

Transport coefficients depend both on the density of states (DOS) and on the scattering 
probability. The behavior of DOS in magnetic field was discussed in Chapter I. We 
have have seen that DOS oscillates because of the energy quantization. The scattering 
probability, in its term, is also dependent on the density of states, as well on the scattering 
matrix element. Consequently, it also oscillates in magnetic field, and it appears that the 
last contribution is the most important. The quantum oscillations of conductivity is called 
the Shubnikov-de Haas effect. Quantum oscillations of transport coefficients are widely 
used for investigation of the properties of metals and semiconductors. 

Let us outline main principles of these effects. To take the electric field into account 
one should analyze the Schrodinger equation in crossed electric and magnetic field (H || 
z, E || x) 

h 2 d 2 i) l fhd e u \ 2 . h 2 d 2 i) , ^ , , n 

H — — — + -Hx ip — -y + (eEx - e)ip = 0. 



2m dx 2 2m \i dy c J 2m dz 2 
Here we use the gauge A = (0, Hx, 0). In this case we can search the solution as 

(p(x) exp(ik y y + ik z z). 

The equation for (p has the form 

h 2 d 2 <p 



2m dx 2 



+ 



— ) x 2 + k y + eE)x + 

2m \ c J \ mc J 2m 



V? = 0. 



The result can be expressed just in the same way as for the case E — with the additional 
terms 

f ffkl , , o m, mc2 ( 1 
e u =e N + — + 5e, 5e — —a H eEk y — I — 

for the energy and 

tp r- 2 r r eEa 2 H 

XQ=xo + bx 0l x = -a H k y , ox = . 

rkd c 

for the oscillator center x$. 

Now we introduce the following concept. Assume that the electron in the state v is is 
situated at the point Xq. The electric current is 

3* = " e E (/o(^) [1 - /o(^)] W% - f (4) [1 - / (ef)] WX] ■ 
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The prime over the sum means that the state v has Xq < 0, while the state v 1 has x^ > 0. 
Then we expand the expression up to the linear in E term and get 

a xx = e^ {—^- ) g W -'- 

v,v' N ' 

This formula has an explicit physical meaning. Indeed, the quantity 

is just the contribution of the states v, v' to the 2D diffusion coefficient in the plane (x, y). 
Thus we come to the old formula 

* = e 2 J deg(e)D(e) (-^f) 

where both g(e) and D(e) should be calculated with the help of quantum mechanics: 

17(e) = D ( £ ) = £ *( e " ^) {X ° ~ X '° r W uv , 

One can see that the result is strongly dependent on the scattering mechanism and oscillates 
in the case of the Fermi statistics. 



2.4 Weak localization 

Consider noninteracting electrons having pp£ ^> H and passing between the points A and 
B through scattering media. The probability is 



W 



5> 



i i^j 



(2.20) 



Here Ai is the propagation amplitude along the path %. The first item - classical probability, 
the second one - interference term. 

For the majority of the trajectories the phase gain, 



Ay? = fir 1 / pdl > 1 , 

J A 



(2.21) 



and interference term vanishes. Special case - trajectories with self-crossings. For these 
parts, if we change the direction of propagation, p — > — p , dl — > — dl, the phase gains are 
the same, and 



1^1 +^| 



\A 1 r+ \A 2 r + 2A 1 A* =4Ui 



2.4. WEAK LOCALIZATION 
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Figure 2.2: Feynman paths responsible for weak localization 



Thus quantum effects double the result. As a result, the total scattering probability at 
the scatterer at the site O increases. As a result, conductance decreases - predecessor of 
localization. 

Probability of self-crossing. The cross-section of the site O is in fact de Broglie electron 

wave length, A ~ fo/pF- Moving diffusively, if covers the distance \J~ij ~ \^Dt, covering the 

volume (Dt) d / 2 b 3 ~ d . Here d is the dimensionality of the system, while b is the "thickness" 
of the sample. To experience interference, the electron must enter the interference volume, 
the probability being 

vX 2 dt 

Thus the relative correction is 




Figure 2.3: On the calculation of the probability of self-crossing. 



Aff r* vX 2 dt 



The upper limit is the so-called phase-breaking time, r^. Physical meaning - it is the 
time during which the electron remains coherent. For example, any inelastic or spin-flip 
scattering gives rise to phase breaking. 



3d case 

In a 3d case, 
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Here we have used the notations 

D~v£, t~£/v, L v = yfD^. (2.24) 
On the other hand, one can rewrite the Drude formula as 

n e e 2 r n e e 2 £ e 2 p 2 F £ 
m pf K 3 

In this way, one obtains 

A "~-s + t (2 - 25) 

Important point: The second item, though small, is of the most interest. Indeed, it 
is temperature-dependent because of inelastic scattering. There are several important 
contributions: 

• Electron-electron scattering: 

T v ~ he F /T 2 . 

Thus 

cx(T)-(x(0) f\\ 3/2 fT\ 



o-(O) \ij \e F J 

It is interesting that at low temperatures this quantum correction can exceed the 
classical temperature-dependent contribution of e-e-scattering which is oc T 2 . It is 
also important to note that the above estimate is obtained for a clean system. It 
should be revised for disordered systems where electron-electron interaction appears 
more important (see later). 

Electron-phonon-interaction. In this case,0 

lip rp3 ' 

and 

a(T)-a(0) f\\ 3/2 fT\ 1/2 f T 



o-(O) \£ J \e F J \huj 



- 



■D 



Thus there is a cross-over in the temperature dependencies. To obtain the dominating 
contribution one has to compare r" 1 . Consequently, at low temperatures e-e-interaction 
is more important than the e-ph one, the crossover temperature being 

T ~ (hw 2 D /ep) ~ 1 K. 



3 Undcr some limitations. 



2.5. WEAK LOCALIZATION .. MAGNETIC FIELD 
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Low-dimensional case 

If the thickness b of the sample is small 0, then the interference volume X 2 v dt has to be 
related to (Dt) d ^ 2 b 3 ~ d . For a film d = 2, while for a quantum wire d = 1. As a result, 



<7 



D 1 b~ 2 L v , d = 1 



It is convenient to introduce conductance as 

G = ab 3 - d . 

We have, 

AG ~ 4 ( ' n < L ^ . (2.27) 

h { L v , d = 1 

Important note: In low- dimensional systems the main mechanism of the phase breaking 
is different. It is so-called low-momentum-transfer electron-electron interaction which we 
do not discuss in detail 



2.5 Weak localization in external magnetic field 

In a magnetic field one has to replace p — ► p + (e/c)A (remember - we denote electronic 
charge as — e). Thus the product AA* acquires an additional phase 

Acp H = ^ /Adl = % /curl AdS = 47rJ- (2.28) 

chj chj $ 

where $ is the magnetic flux through the trajectory while $o = 2nhc/e is the so-called 
magnetic flux quantum. Note that it is 2 times greater than the quantity used in the theory 
of superconductivity for the flux quantum. 

Thus magnetic field behaves as an additional dephasing mechanism, it "switches off" 
the localization correction, increases the conductance. In other words, we observe negative 
magnetoresistance which is very unusual. 

To make estimates, introduce the typical dephasing time, to get Aipn ~ 27r for the 
trajectory with L ~ y/Dtfj. In this way, 

t H ~ %/(HD) ~ l 2 H /D , l H = y/ch/eH . (2.29) 

The role of magnetic field is important at 



4 The general from of the critetion depends on the relationship between the film thickness, b, and the 
elastic mean free path £. The result presented is correct at b <C I. At L 3> b^$> £ one can replace the lower 



limit t of the integral ( p.22Q by r fc = 6 2 /D 
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Note that at H « H 

u c t ~ h/(e F T v ) < 1 

that means the absence of classical magnetoresistance. Quantum effects manifest them- 
selves in extremely weak magnetic fields. 

To get quantitative estimates one has to think more carefully about the geometry of 
diffusive walks. Let consider the channel of 2DEG with width W . The estimates given 



above are valid for £, L v <C W, exact formulas could be found, e. g. in Ref. [10]. Usually 



Hq is very weak, at L v — 1 /im Hq w 1 mT. The suppression of weak localization is 
complete at H > h/e£ 2 , still under conditions of classically weak magnetic field, uj c t <C 1. 

The situation is a bit different at W < this case can be mentioned as one- 
dimensional for the problem in question, see Fig. |2.4| . Then a typical enclosed area is 




w 



Figure 2.4: On quasi ID weak localization. 



S ~ WLp, and the unit phase shift takes place at 



t H ~ L%/DW 2 , 



Some experimental results on magnetoresistance of wide and narrow channels are shown 
in Fig. ^75| . 

There are also several specific effects in the interference corrections: 

• anisotropy of the effect with respect of the direction of magnetic field (in low- 
dimensional cases); 

• spin-flip scattering acts dephasing time; 

• oscillations of the longitudinal conductance of a hollow cylinder as a function of 
magnetic flux. The reason - typical size of the trajectories. 



2.6 Aharonov-Bohm effect 



Magnetoconductance corrections are usually aperiodic in magnetic field because the in- 
terfering paths includes a continuous range of magnetic flux values. A ring geometry, 
in contrast, encloses a continuous well-defined flux $ and thus imposes a fundamental 
periodicity, 

= G($ + n$ ) , % = 2nhc/e. 



2.6. AHARONOV-BOHM EFFECT 



29 




8(G) 

Figure 2.5: Experimental results on magnetoresistance due to 2D weak localization (upper 
panel) and due to ID weak localization in a narrow channel (lower panel) at different 
temperatures. Solid curves are fits based on theoretical results. From K. K. Choi et ai, 
Phys. Rev. B. 36, 7751 (1987). 



Such a periodicity is a consequence of gauge invariance, as in the originally thought exper- 
iment by Aharonov and Bohm (1959). The fundamental periodicity 

eS 

comes from interference of the trajectories which make one half revolution along the ring, 



see Fig. |2.6| . There is an important difference between hc/e and hc/2e oscillations. The 
first ones are sample-dependent and have random phases. So if the sample has many 
rings in series or in parallel, then the effect is mostly averaged out. Contrary, the second 
oscillations originate from time-reversed trajectories. The proper contribution leads to a 
minimum conductance at H = 0, thus the oscillations have the same phase. This is why 
/ic/2e-oscillations survive in long hollow cylinders. Their origin is a periodic modulation of 
the weak localization effect due to coherent backscattering. Aharonov-Bohm oscillations 
in long hollow cylinders, Fig. [2.7| , were predicted by Altshuler, Aronov and Spivak [EO] and 
observed experimentally by Sharvin and Sharvin [21]. A rather simple estimate for the 
magnitude of those oscillations can be found in the paper by Khmelnitskii [22]. 
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Figure 2.6: Illustration of Aharonov-Bohm effect in a ring geometry, (a) Trajectories 
responsible for hc/e periodicity, (b) trajectories of the pair of time-reversed states leading 
to /ic/2e-periodicity. 




Figure 2.7: On the Aharonov-Bohm oscillations in a long hollow cylinder. 



2.7 Electron-electron interaction in a weakly disor- 
dered regime 



Let us discuss the effect of the e — e interaction on the density of states. Let us concentrate 
on the exchange interaction, shown in a right panel of Fig. ^T5| 



Ae 



<7(k) 



|p-frk|<p F 



d 3 k 



(2.30) 



2.7. E-E INTERACTION 
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Here g(k) is the Fourier component of the interaction potential, the sign "-" is due to 
exchange character of the interaction. In the absence of screening g(k) = 4ne 2 /k 2 , and 



Ae 



kr 



d 3 k 



|p-?ik|<p F 



-4vre 2 /i 2 

„2 



P'<PF 
2 



(2tt) 3 
1 d 3 p' 
|p -p'| 2 (2nhY 



Pf 



p-pp P + Pf 



2p 



In 



P-PF 



To obtain this formula it is convenient to use spherical coordinates - (dp') 
and auxiliary integral 



d(cos 0) 



i p2 _|_ p/2 _ 2pp' COS ( 



1 

ppl 



In 



p + p' 
p — p> 



(2.31) 



2np' 2 dp' d(cos 0) 



(2.32) 



At small p — pf it is convenient to introduce £ = v{p — pp) -C ep to get (omitting the 
irrelevant constant) 

Ae = -(e 2 /nhv)£ \n(2p F v/£) . (2.33) 

Screening can be allowed for by the re place ment k~ 2 — > (k 2 + n 2 )" 1 that replaces p±p' in 
the argument of the logarithm in Eq. (2.32) by \J(p ± p') 2 + (hn) 2 . As a result, 

-(e 2 /nhv)^\n(2pp/hK) at Hk^Pf- 



Ae 



(2.34) 



This is a simple-minded estimate because it ignores the interference of the states. Indeed, 
if the two states differ in the energy by |£| the coherence time is h/\£\. If the electron 
returns back for this time, then the effective interaction constant increases by 



Thus 

In a similar way 



C eff 

Av _ 

v 



v\ 2 dt 



e 2 (l + a € ). 



d ' 



(2.35) 



Ag 

9 



hv 
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Here g is the density of electronic states. Consequently, since a oc v, we arrive to a specific 
correction to conductance. Comparing this correction with the interference one we conclude 
that interaction dominates in 3d case, has the same order in 2d case and not important 
in Id case. Another important feature is that the interaction effects are limited by the 
coherence time h/\£\ ~ r T = h/k B T rather than by r^. Usually, r v 3> t t . Consequently, 
the interference effects can be destroyed by weaker magnetic fields than the interaction 
ones (important for separation of the effects). 

2.8 Few words about Anderson localization 

If we come back to the interference correction for d — 2, 1 we observe that it increases with 
t v , or at T — > 0. Thus the corrections becomes not small. We can also prepare the samples 
with different values of the mean free path i. 

What happens if the corrections are not small? Anderson (1958) suggested localization 
of electronic states at T = 0. This suggestion has been later proved for an infinite Id 
system, as well as for an infinite wire of finite thickness (Thouless, 1977). Later it has been 
shown that 



where L loc ~ i in the first case and (bpp/h) 2 £ for the second one (exponential localization). 
It seems that such a law is also the case for a metallic film (rigorous proof is absent). 
Very simple-minded explanation - over-barrier reflection + interference of incoming 
and reflected waves. Because of the interference the condition T = is crucial (no phase 
breaking). This explanation is good for one-dimensional case. 

A little bit more scientific discussion. Consider interference corrections to the con- 
ductance at T — > 0. In this case one has to replace 



One can conclude that in 3d case the relative correction is ~ (H/pf^) 2 <C 1 (usually). 
However, at d — 1, 2 it increases with the size. 
At what size A<r ps a? 



We can conclude that in 3d case localization takes place at pp£ ~ H, while in Id and 2d 
case it takes place at any concentration of impurities. 

Scaling hypothesis - According to the "gang of 4" (Abrahams, Anderson, Licciardello 
and Ramakrishnan) 



G oc exp(-L/L loc ) 



L 2 /D, L^^L. 




(2.36) 





G(L) = f[l,G(L)}; 
aLG'(L) = a(df/dq) q=1 . 



(2.38) 



2.8. ANDERSON LOCALIZATION 
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Denoting 

G-\df/dq) q=l = (3(G) 
we re- write the scaling assumption through the Gell-Mann & Low function, (3(G): 

d\nG/dlnL = (3(G). 

At very large G we can expect that the usual theory is valid: 

( Sl/L = L, d = 3; 
G = a\ bL ± /L = b, d = 2- (2.39) 
{ b 2 /L, d = 1 

Thus, in the zero approximation, 

(3(G) = d-2. 

Then we can use the weak localization approximation to find next corrections. One can 
get 

(3(G) « d - 2 - a d G /G , (2.40) 

where 

G = e 2 / (ix 2 K) , a d ~ 1 . 

Indeed, for d = 3 

In G = \n[(a + Aa)L] ^\naL + (Aa)/a w In [(a + Act| l=00 ) L] + h 2 /(p 2 F £L) . 

Thus 

<91nG ^ 2 n 2 a G 

m= d^ = 1 'p^L = 1 -p 2 UG = 1 - a3 ^- 

At small G one can suggest exponential localization: 



G ~ G exp(-L/L c ) -> /3(G) ~ ln(G/G ) . 

Thus we have the scenario shown in Fig. ^!.9| . Believing in such a scenario we get localization 
for Id case ((3(G) < - conductance increases with the length). At d = 3 we have a fixed 
point at G c , which is unstable ((3 changes sign). Under the simplest assumptions 

(3(G) « t (1iiG-1iiG c )«7(G/G c -1), 
G = G 1 - ^ at L = Lq (initial condition) 

(G^ is close to G c ) we obtain 
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From Eq. (|2.41 ) one can find important dependencies near the critical point. It is natural 



to chose Lq ?a I and to suggest that at this size 

~ e 2 p 2 F e/h 3 -> G = a L = (e 2 /h) (p F £/h) 2 . 

Now we can assume that we can control some parameter (say, impurity content, x, which 
effects the mean free path £), and that Go is regular in this parameter. Then, at small x 
= G c (l + x). At L = i we obtain 

G = G c (l + « G c exp[x{L/ey] . 

Of course, such an assumption valid only at 

x < 1 , LI I > 1 . 

Thus at x < we obtain exponential localization with the characteristic length 

L loc ~£|x|- 1/7 , x<0 - 

However, at x > G grows with L. In a spirit of the concept of phase transitions we can 
treat the quantity L c = tx~ x ^ as a correlation length. At the scales of the order of L c the 
properties of conducting and insulator phases are similar. The above law can be valid only 
in the vicinity of Go ~ G c . Then we match the usual Ohm's law, G = o~L. Consequently, 
the conductivity can be estimated as 

Go / G \ x i„ 

Thus, we predict power law. Experiments on the dielectric constant (re oc L 2 ) support the 
value 

7 = 0.6 ±0.1. 



36 



CHAPTER 2. DIFFUSIVE TRANSPORT 



Don't forget that we discuss the case T = 0, and the conductance is supposed at zero 
frequency. The range of applicability is given by the inequalities 

L c < L v , = VTJuj . 

This is almost impossible to meet these conditions, so usually people extrapolate experi- 
mental curves to T = 0, u — 0. It is a very subtle point because, as it was shown, the 
conductance at L < is not a self- averaging quantity with respect to an ensemble of 
samples. More precise, the fluctuations between the samples cannot be described by the 
Gaussian law, their distribution being much wider. Then, 

• have all these scaling assumptions any sense? 

• Why they reasonably agree with the experiment? 

The answer is positive because both the scaling predictions and the experiment are valid 
as an extrapolation from the region L > L v . 

As the temperature grows, fluctuations decrease and the conductance tends to the 
Ohm's law. 

Role of e — e interaction 

Nobody can consider both disorder and interaction acting together. To get some under- 
standing of the role of e — e interaction let us consider a clean metallic conductor. Assume 
that under some external perturbation (say, pressure) the band overlap changes. In this 
way we control the Fermi level (number of electrons and holes). One can consider them 
as free ones as their kinetic energy p 2 F /2m* exceeds the potential energy e 2 /(/t f). In this 
way we come to the condition 



KqHv 

Otherwise electrons and hole form complexes - Wannier-Mott excitons. This state is insu- 
lating because excitons are neutral. This is only one of possible scenario. In general, the 
problem is a front end of modern condensed matter physics. 



Chapter 3 
Ballistic transport 



3.1 Landauer formula 

We start this chapter by a description of a very powerful method in physics of small systems 
- so-called Landauer approach. 

The main principle of this approach is the assumption that the system in question is 
coupled to large reservoirs where all inelastic processes take place. Consequently, the trans- 
port through the systems can be formulated as a quantum mechanical scattering problem. 
Thus one can reduce the non-equilibrium transport problem to a quantum mechanical one. 

Another important assumption is that the system is connected to reservoirs by ideal 
quantum wires which behave as waveguides for the electron waves. We start our analysis 
from the discussion of the properties of an ideal quantum wire. 



Ideal quantum wire 

Consider 2 large reservoirs of electron gas reservoirs having the difference Sn in the electron 
density and separated by a pure narrow channel. For small 8n one can assume that there 
is a difference in a chemical potential, 5/i = 5n/g{ep)- In the following we shall use the 
Fermi level of non-biased system as the origin for the chemical potentials. So the difference 
between the chemical potential in ct-th reservoir will be denoted as fi a . 

If the channel is long and uniform, then the total current carried by the state char- 
acterized by a transverse mode n and a given direction of spin which propagates without 
scattering is 

j =e r dk z de n (k z ) = 2 r F+ ^ dc de n (k z )/dk z = 2 5 
J 2nh dk z 2nh J ep+fla " \de n {k z ) / dk z \ h 

If we take into account electron spin and N transverse modes are open, then the conduc- 
tance is given by the expression G = H-N. 

We come to a very important conclusion: an ideal quantum wire has finite resistance 
h/2e 2 N which is independent of the length of the wire. 
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As we have seen, even an ideal quantum wire has a finite resistance. That means a 
finite heat generation even in the absence of any inelastic processes inside the wire. Below 
we will discuss the physical picture of heat release by a current-carrying nanostructure 
(here we follow the considerations of Ref. [E3]). 

First of all let us specify what heat release is. It will be convenient to consider an 
isolated system. Therefore we will have in mind the following physical situation. There is 
a capacitor which is discharged through the conductor of interest. The product RC of the 
whole system, R and C being the resistance and capacitance respectively, is much bigger 
than any relaxation time characterizing the electron or phonon system of the conductor. 
This means that for all the practical purposes the conduction process can be looked upon 
as a stationary one. The total energy of the system, U, is conserved, while its total entropy, 
S, is growing. The rate of heat generation is expressed through TdS/dt, where T is the 
temperature, i.e. through the applied voltage and characteristics of the nanostructure itself. 
This means that the result is independent of the assumption that the considered system 
is isolated, which is made only for the sake of derivation. This thermodynamically defined 
heat is generated in the classical reservoirs over the length having a physical meaning of 
the electron mean free path. That is the same mean free path that enters the Drude 
formula, which determines the conductivity of the reservoirs themselves, the amount of 
heat generated per second in both reservoirs being the same. 

It is interesting to indicate that even purely elastic collisions can result in a heat gen- 
eration although they of course cannot establish full equilibrium. This has a clear physical 
meaning. The amount of order in the electron distribution resulting in electric current can 
bring about mechanical work. For instance, one can let the current flow through a coil, 
and a magnetic rod can be drawn into the coil. In such a way the electrons transferring 
the current can execute a work on the rod. As a result of scattering, the amount of order 
in the electrons' distribution diminishes, and this means dissipation of mechanical energy 
into the heat. It has been shown that the heat release is symmetric in both reservoirs even 
if the scatterers in the system are asymmetric. 

All the above considerations do not mean that the collisions that give the main con- 
tribution to the heat release, also establish full equilibrium. What equilibrium needs is 
inelastic collisions which transfer the energy of electrons taking part in charge transfer to 
other degrees of freedom, such as to other electrons and phonons. In particular, a local 
equilibrium electron distribution is established over the length scale determined by electron- 
electron interaction. Such a distribution can be characterized by a local electro-chemical 
potential and sometimes an electron temperature. The latter can in principle be mea- 
sured by optical methods. On the other hand, the equilibrium with respect to the lattice 
is established at the scales of electron-phonon and phonon-phonon mean free paths. Only 
over those distances from the channel one can treat the results in terms of the true local 
temperature. 
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Resistance of a quantum resistor 

Consider a system shown in Fig. |3.1| consisting of a barrier connected to reservoirs by ideal 
quantum wires. If there is some reflection only a part of the current is transmitted. In this 




Figure 3.1: On the resistance of a quantum resistor. 

case one can introduce the transmission probability of the mode n, T n , to obtain (including 
spin degeneracy) 

2 N 
J= h 5fl ^ Tn ' 

n=l 

As a result, 

9 P 2 JL op 2 
G=^£T n = ^Trttt. (3.1) 

n=l 

Here t is the matrix of scattering amplitudes while the expression is called two-terminal 
Landauer formula. 

This very important and looking simple formula was confusing during a long period. 
Indeed, this is the conductance which is measured between two reservoirs. Having in mind 
that the resistance of the connecting ideal wires (per one conducting mode) is h/2e 2 we 
can ascribe to the scattering region the resistance 



h 

2^2 



1 

r" 1 



h R 
2^2 f 



where R is the reflection coefficient. Consequently, in the original formulation the quantum 
resistance was described as 

2e_ 2 A T n 
h 



a = - E rzV • (3 - 2) 



n=l 

However, the quantity which is usually measured is given by Eq. (|3.1| ). 

Now we derive the Landauer formula for finite-temperature and so-called multichannel 
case when the leads have several transverse modes. Consider ideal wires which lead to 
a general elastic scattering system. Let each lead has the cross section A and have N± 
transverse channels characterized by wave vectors ki so that, 

h 2 k 2 
2m 
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The incoming channels are fed from the electron baths with the same temperature and 
chemical potentials f/,i, f/,2, ■ ■ ■■ The outgoing channels are fed up to thermal equilibrium 
population. We shall assume that the particles are absorbed in the outgoing baths. The 
sources are assumed to be incoherent, the differences — /i 2 are also assume small to yield 
linear transport. We introduce the scattering amplitudes for the transmission from jth 
incoming to ith outgoing channel. Reflection amplitudes are introduces in a similar 
way for reflection into the ith incoming channel. If we replace the incoming and outgoing 
channels, we denote the proper amplitudes by primes. In this way it is convenient to 
introduce 2N± x 2N± scattering matrix as 

r t' 
t r' 

From the current conservation we must require unitarity while from time reversal symmetry 
S — S. Thus we have also SS* = I where star stays for complex conjugation while tilde 
for transposition. In a magnetic field the Onsager relation requires S(H) = S(—H). 
It one defines the total transmission and reflection into ith channel as 

Ti = ^ ] \tij\ , Ri = ^ ] \ r ij\ ■ 
j j 
then from unitarity condition we get 

i i 

Since the densities of states in each channel are ID like, gi(E) = (nhvi)" 1 we write the 
current through outgoing channels as 



e 

7ih~ 



Thus the conductance becomes 



Y.J dE W E W E ) + h(E)R'(E) - ME)] 

( " i: " 2) 7 d£ (-i)l> (£) ' 



h J \ dE 

This is the two-terminal conductance measured between the outside reservoirs which in- 
cludes contact resistances. 



Multiterminal resistance 

For simplicity we shall discuss the case of zero temperature. Let us introduce the total 
transmission probability from the bath a to the bath f3, 

T a — >p ^ ^ ^ ^ |^/3a,mn| • 
n=l m=l 
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Here N{ is the number of propagating modes in each lead connected to zth reservoir. 
Counting all the chemical potentials from the Fermi level, we see that the reservoir a 
injects the current {2ejh)N a \i a into the lead a. The fraction T a _>p/N a is transmitted to 
the reservoir j3 while the fraction T a ^ a /N a = R a /N a is reflected back into reservoir a. 
The the net current I a is given by the following set of equation, 



h 
2^ 



(3-3) 



Introducing vectors I and jl with components I a and fi a , respectively, we can write 

1 = 6(1, (3.4) 

where the conductance matrix G is defined as 

2e 2 



G 



ai3 



h 

h 



[(N a - R a ) 5 af3 - Tp_+a (1 - 6 a/ 3)] 
[N a 5 a p — Tp_> a ] . 



(3.5) 



Here we use the relation T a ^ a = R a . The sum of rows of this matrix is zero because 
of current conservation, the sum of the elements of each row also vanishes because if one 
changes all the chemical potentials by the same amount no current will be induced. Thus 

N a — R a = Tp^ a = T a ^/3 . 

are called often the Landauer-Biittiker formalism. They 



The equations 



and 




Figure 3.2: On the resistance of 4-terminal device. 

allow find, e. g. 4-terminal resistance. Indeed, we can put I\ = —J; 
Then / = Ij where 

/ 1 \ 

-1 

J 



/, h 



Thus 



n 
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Having in mind the properties of the scattering amplitudes we have, 




W-#) 



that results in the reciprocity relation 



(H) 



Here lZ a p tl s stands for the resistance measured for voltage contacts 7, 5 while the current 
passes through the contacts a, j3. Note that this relation works even in the case when the 
concept of local conductivity is not applicable. What we only need is linear response and 
absence of inelastic scattering inside the device under consideration. 

One can easily generalize the above expressions for the case of finite temperatures by 
replacement of the element of G-matrix by their thermal averages, 



3.2 Application of Landauer formula 
Point ballistic contact 

The most clean system is the so-called quantum point contact (QPC) - short and narrow 
constrictions in 2d electron gas. A sketch of QPC is shown in Fig. |0| The conductance of 



QPC is quantized in the units of 2e 2 / h. The quantization is not that accurate as in quantum 
Hall effect (about 1%) because of non-unit transparencies T n and finite temperature. It is 
interesting to compare quantum and classical behavior of QPC. In a classical picture one 
can write 



(A) T 



J °°dEg(E) (df /dE)A(E) 
j °°dEg(E) (df /dE) 




Figure 3.3: A sketch of QPC formed by splitted gates. 




Thus the "diffusion constant" is 



D, 



J _ 1 

on 7r 



Wv F 



G = e 2 g(e F )D eff 



2e 2 k F W 
h 71 
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Gate voltage V a [V] 



Figure 3.4: Quantization of conductance of a point contact: Schematic picture (left) and 
experimental result (right). 



D 




Figure 3.5: On the classical conductance of a point contact. 

Note that the integer part of the quantity k F W/n is just the number of occupied modes 
according to quantum mechanics. 

Series addition of quantum resistors 

Assume that we have two obstacles in series. Let the wave with unit amplitude is incident 
to the region, the amplitude of the reflected wave is A while D is the amplitude the 
wave transmitted through device. The obstacles are connected by an ideal conductor, the 
phase shift of the wave along which being 0. Let the wave emerging from the obstacle 



1 




id) 

B Be Y 




D 








A 


C Ce^* 



Figure 3.6: On series of quantum resistors. 



1 is B exp(kx — out). It reaches the obstacle 2 gaining the phase 0, having the complex 
amplitude Bexp(i(p). The reverse wave C gains the phase —0. In this way we get the 
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following set of equations, 

A = n + t]C , B = t! + r[C 

Ce~^ = r 2 Be i4, , D = t 2 Be ilt) 

Solving this equation we obtain, 

D 



1 - e 2i <t>r 2 r' 1 ' 
that yields for the total transmittance: 

T T 

T=\D\ 2 = 1 % (3.6) 

1 + R 1 R 2 -2 v ^R 2 cose 

where 9 = 20 + arg(r 2 r' 1 ). The ratio between reflection and transmission which should be 
understood as a reduced resistance (in units h/2e 2 ) of the system excluding wires is 



n _ x _R A 2 R l + R 2 -2^Rjl 2 cos6 „ _ KG 
G =f = 



T{T 2 ' 2e 



G = (3-7) 



This is a very strange formula. Assume that we made an ensemble of the systems which 
differ only by the distance between the obstacles, i. e. by the phase 0. Let the distribution 
of (j) will be constant in the interval (0, 2ir). The averaging over we get 



[1 - R 1 )(l - R 2 ) ' 
while the Ohm's law will provide 

Ri R 2 



Q- l = 



1 - R-y l-R 2 



As a result the Ohm's law survives only at small reflections. 

Let us construct a chain of n resistors with very small reflections. Then the total 
reflection first increases linearly in n. Finally the total transmission becomes substantially 
less than 1. Now let us add a very good conductor to this chain. We get 

{y )n+i = — = {y )n + 7=r- 

Thus an addition a good conductor increases the resistance by R/T n > R. Such a behavior 
can be formulated as a "renormalization group" 

s>->. = «r'>„ + i. 

Thus the average resistance grows exponentially with the length which has something to 
do with ID localization. This considerations are not fully satisfactory because resistance 
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is not the proper quantity to be averaged. Following Anderson, the proper quantity to be 
averaged is ln(l + Indeed, 



1 + = 1 + R/T = 1/T ^ln(l + £- 1 ; 



InT. 



The quantity — In T plays the role of extinction exponent and it should be additive to 
successive scatterers is the relative phases are averaged out. We get this relation using 

r 27T i r 

/ dO ln(a + bcosO) = Trln - a + -y/a 2 — b 2 
Jo 2 . . 

So the exact scaling is given by the relation 

(ln(l + g- 1 ))=n(Rh/2e 2 ). 

Parallel addition of quantum resistors. 

Let us now discuss the parallel addition of two single- channeled quantum resistors. The 
geometry of the problem is shown in Fig. |3.7| . All the phases and scattering effects along 
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Figure 3.7: On the parallel addition of quantum resistances. 

the branches are absorbed by the scattering parameters. Time-reversal symmetry requires 
ti = t' { while the current conservation requires 

-U/t't = n/r't . 

In the presence of Aharonov-Bohm flux $ through the loop, following from the gauge 
invariance the scattering amplitudes are renormalized as 

h -> t ie - * t[ -> t ie +ie , t 2 -> t 2 e +i6 ,t' 2 -> t 2 e- ie , n -> r h r: 



r' 



Here 6 = 7r$/$ - 

This point needs some more explanation. An Aharonov-Bohm flux $ through the 
opening can be represented as § A • dl along the path circulating the opening. Here A is 
the vector potential. One can eliminate this flux by a gauge transform 



if)' = if) exp 



he 
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where \ is defines as Ai = V%. The transformed Schrodinger equation has A\ = 0. 
However, the price for this is that the transformed wave function, ip', does not satisfy 
periodic boundary conditions. When the electron coordinate is rotated once around the 
ring the phase of x' is changed by 5\ = 27r$/$ - in our calculation this phase shift is 
absorbed into the expressions for the transition amplitudes. 

To find the transmitted wave one has to determine 10 unknown amplitudes, Xi,x 2 ,yi,y 2 ,Ui,u 2 ,Vx,v 2 , F, 
It can be done solving the set of matching equations at the scatterers and the triple con- 
nections. It is assumed that the connections do not introduce additional scattering and 
can be described by the unitary scattering matrix 

/ -1/V2 -1/V2\ 

S = -1/V2 1/2 -1/2 . 
V -1/V2 -1/2 1/2 / 

Here Sa denote the reflection amplitude of the ith channel while off-diagonal elements Sij 
are the transition amplitudes from the channel i to j. The subscript 1 is used for left 
incoming channel and right outgoing channel. After rather long algebra made originally in 
Ref. [E33] we arrive at the solution, 

T-\F\ 2 = i a + f 3cos26 ( oo) 
1 1 7 + 5cos2# + ecos4# ' { ' } 

where a, j3, 7, 5, e are rather complicated functions of the scattering amplitudes, 

a = \A\ 2 + \B\\ f3 = 2U(AB*), 7 = \D\ 2 + \E\ 2 , 
5 = 2^ (DC* + EC*) , e = 2 3? (DE*) , 

A = tjt 2 + t 2 (n - 1)(1 - r[) , B = tjh + tx(r 2 - 1)(1 - r' 2 ) , 
D = E = ttt 2 C = t\ + 1\ - (2 - r x - r 2 )(2 - r[ - r' 2 ) . 

This expression describes a rich physical picture. Even in the absence of magnetic field, 
6 = the transmittance can be strongly dependent on the phases of the complex scattering 
amplitude. If we make one branch fully non-conducting, t\ = 0, still the appropriate choice 
of phases of the reflection amplitudes r± and r[ can result either in T = or T = 1. Indeed, 
in this case, 

A = t 2 {n - 1)(1 - ri) , B = D = E = , C = t\ - (2 - n - r 2 )(2 - r[ - r' 2 ) , 
a=|t 2 | 2 |(r 1 -l)(l-ri)| 2 , P = 0, 7 = \t\ - (2 - n - r 2 )(2 - r\ - r' 2 )\ 2 , 
5 = e = 0. 

As a result we get 

T 



t 2 \ 2 \(n - - r[) 



|^-(2-r 1 -r 2 )(2-r / 1 -r 2 )| 2 ' 

Putting ri = 1 we obtain T = 0. Thus we observe that non-conducting branch can influence 
the total conductance strongly. 

Of course, all the discussed effects are due to interference. If the size of the system 



exceeds L v we come back to classical laws. 



3.3. ADDITIONAL ASPECTS OF BALLISTIC TRANSPORT 



47 



3.3 Additional aspects of ballistic transport 
Adiabatic point contacts 

The results of first observations of conductance quantization were surprising. Indeed, from 
quantum mechanics it is well know that any sharp potential barrier produces oscillations 
in the transmission coefficient as a function of the energy. However, the experimental steps 
were more or less rectangular. An explanation of such a behavior was given in Ref. [E3]. 
The authors showed that if the point contact has a smooth profile, i. e. if its width d 
depends on the longitudinal coordinate x in a smooth way, then the T(E) dependence is 
very close to a perfect step. To make the results simple, let us consider a channel with 
rectangular confinement. Let us assume that we can separate the variables in an adiabatic 
way, ty n (x,y) = ip n (x)(p njX (y), and first solve the Schrodinger equation for a given width 
d. In this way we get the transverse wave functions 



<fn,x{v) 



d(x 



sin 



7m- 



2y + d(x) 



d(x) 



The Schrodinger equation for the longitudinal motion has the form 



k 2 d 2 i) 
2m dx 2 



+ e n (x)ip = Eip , e n (x) 



7T 



2 n 2 k 2 



2m[d(x)] 2 ' 

If the variation d(x) is smooth at the scale of de Broglie wave length, kp 1 , the potential 
e n (x) is semiclassical. Then one can use the semiclassical scheme for scattering problem 
and choose 



IPnXOO) 
Pn{x) 



exp 



h 



Pn(x') dx' 



Pn{X) 



'2m[E - e n {x)\ . 



The transmittance step depends occurs when the Fermi energy crosses the maximum of 




| Energy 



Figure 3.8: On the adiabatic quantum point contact. 



the potential e n (x) for the upper transverse mode, see Fig. |3~8| . Expanding the potential 
near its maximum we get 

. . [ / d 2 d(x)\ x 2 ' 
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Since d 2 d/dx 2 — 2/p where p is the curvature radius of the center of constriction, we get 
the barrier as 

T 2x 2 ~ 
U(x) = e n (0) l-^- 
dp 

The transmission through a parabolic barrier is known, 

1 



T{E) 



1 + exp[-7r 2 (A;do/7r - n )y / 2p/d \ 



(3.9) 



Here do is the minimal width of the constriction, uq is the number of upper level, while 
k = hr 1 \/2mE. We observe that the shape of the step is almost n independent, the 
transition being sharp at p > d . It is important that the numerical factor n 2 \/2 makes 
the transitions sharp even at R ~ d . The same numerical factor helps for the semiclassical 
condition to be valid. This criterion reads n 2 y^2p/d( ) 3> 1. To make the motion through 
the contact ballistic the elastic mean free path should exceed y/ pd . 



3.4 Electron-electron interaction in ballistic systems 

The case of pure 3D metal. Concept of Fermi liquid 

Let us begin with the estimate of the electron-electron scattering in a Fermi gas. Suppose 
that we have a particle 1 outside the Fermi sea, see Fig. EH| If this particle interacts with 




Figure 3.9: Scattering processes for electron-electron interaction. 

another one, 2, inside the Fermi sea both final states should be outside the Fermi sea (Pauli 
principle!). According to the momentum conservation law, 

Pi + P2 = Pi + P' 2 , 

and, as we have seen, 

Pi, Pi, P2 > Pf] P2<Pf. 

The momentum conservation law is shown graphically in the right panel of Fig. |3.9| . One 
should keep in mind that the planes (pi,P2) and (Pi,p' 2 ) are not the same, they are 
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shown together for convenience. To get the escape probability for the particle 1 one should 
integrate over the intermediate momenta 



W oc J 5{ei + e 2 - e\ - e' 2 ) (dp 2 ) (dpi) 



(p' 2 is fixed by the momentum conservation). The energy conservation law actually deter- 
mines the angle between p' : and p' 2 for given absolute values of these vectors. Consequently, 
the rest is to integrate over p 2 = \p 2 \ and p[ = |px|. 

Let pi be close to pp. It means that all the momenta are close to p F and the angles 
with the vector pi + p 2 are almost the same. So let us assume cosines to be the same and 
from the relation between the projections write down 

p\ ^pi+p 2 -p' 2 - 

Now let us recall that p' 2 > pp. Consequently, p[ < pi + p 2 — pp. But at the same time, 
Pi > Pf- Thus 

Pi + Pi ~ Pf > Pf, or p 2 > 2p F - p 1 . 

But from the Pauli principle, the upper limit for p 2 is p F . As a result, we come to the 
following chain of inequalities 

> p 2 - Pf > Pf - Pi, < p\ - p F < (p ± - p F ) + (p 2 - p F ). 

Finally, 

/f° rai+a-2 ^2 

dp 2 dp\ = / da 2 / da[ = — - 
J- ai Jo 2 

where we have introduced cti = Pi—p F . Now we should remember that e — e F = v F {p — p F ). 
So W oc (e— e F ) 2 . The simplest way to estimate r is to use dimensionality approach. Indeed, 
the average potential and kinetic energies are of the order of e F . Consequently, the only 
quantity which is proportional to (e — e F ) 2 and has the time dimensionality is 

he F 

T ~ 



{e-epf 

We came to important conclusion: if one is interested in quasiparticles near the Fermi 
level, \e — e F \ <^ e F , he can treat them as to near classical ones provided 

h e — e F 

< 1. 



[e-e F )T e F 

The typical value for the quasiparticle energy is k^T. This is why the electron-electron 
interaction can be treated in the leading approximation in a self-consistent approximation. 
The estimates above based on the conservation laws provide background of the theory 
of Fermi liquid. The essence of this concept is that the excitations in the vicinity of 
the the Fermi surface can be treated as quasiparticles which behave as particles with 
renormalized velocity. Consequently, the effects of electron-electron interaction are not 
crucially important in pure 3D systems. 
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One-dimensional systems. Tomonaga-Luttinger liquid 

For ID interacting systems the above considerations are not valid because for a single 
branch linear dispersion near the Fermi points the energy spectrum is close to linear, 
E — Ep pa v(p — pf). That means that the energy and momentum conservation laws 
are actually the same, and this is why they are not restrictive as in a 3D case. For this 
reason, the perturbative corrections describing even weak electron-electron interaction are 
divergent. A proper model for interactive ID electrons is the so-called Tomonaga-Luttinger 
model. According to this model, collective electron modes (plasmons) with linear spectra 
are described by new, boson modes. Creation of a real electron in this model is equivalent to 
excitation an infinite number of plasmons. Because of that, the space and time dependence 
of density (and spin) correlation functions are substantially different from the ones for non- 
interacting systems. That manifests itself in various kinetic quantities. For example, the 
Drude conductivity is predicted to vary as power law with temperature. 

The Luttinger liquid model which was previously used for ID organic conductors now 
became important for high-mobility quantum wires, as well as for edge states under con- 
ditions of quantum Hall effect (see below). One can find a good review of this model in 



A ID quantum wire is appropriately characterized by a conductance. It the absence 
of interactions, the conductance of an ideal single-mode quantum wire, adiabatically con- 
nected to leads, is quantized, G = 2e 2 /h. In the presence of a scatterer, the conductance 
drops to G — 2e 2 T/h, where T is the transmission coefficient. 

The electron-electron interaction modifies dramatically the low-energy excitations in 
a quantum wire that leads to striking predictions for the transport. The new features 
manifest itself only if there is one (or several) scatterers inside the quantum wire - otherwise 
the correlation effects are canceled out at the contacts between the interacting quantum 
wire and non-interacting reservoirs. All that together leads to a rich and very interesting 
physical picture. 

To get a flavor of the theory [] let us consider a spinless electrons hopping on ID lattice 
with the Hamiltonian 



When the interaction V = this Hamiltonian can be diagonalized as E^ = —tcosk, \k < 
7r. The low-energy excitation exist near ±kp. Consider a single particle excitation near +kp 
where we remove one electron with k < kp and place it into a free state with k + q > kp. 
Then the energy of excitation is foujk = hqvp. Adding a similar state near — kp we have 
a situation similar to phonons in one dimension. When the interaction is turned on this 
dispersion law remains, however the velocity is renormalized. 

Linear spectrum implies a boson-like description. Mathematically in can be done using 

1 Here we follow Ref. JIB]. 



Ref. [[El]. 




(3.10) 
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(Jordan-Wigner) canonical transform, 



exp | in 

k>j 



C W ) bj , 
k>j J 



that keeps the Hamiltonian (|3.10|) in the same from with replacement c —>■ b. One can 
check that the b operators at different lattice points commute, and therefore they are bosons. 
Now, the boson operators can be approximately decoupled as bj — > v /rTJexp(i0 ? -) , rij = CjCj 



^3 

Then we can proceed to continuum limit, focusing on scales long compared to the lattice 
constant. In this way we replace 

<Pj -> <i>fa) , rij -> p(x) . 

Extracting from the total electron density its average value, po = kp/n and introducing the 
"displacement" operator, 9, as p — p = d x 9(x)/ir we arrive at the phonon-like commutation 
rule, 

ITT 

[cf>(x),6(x')} = ^5(x-x'). 

We observe that d x <j) is the momentum conjugate to 9. As a result, we arrive at the effective 
Hamiltonian, 

n = ^[g(d x <P) 2 + g-\d x 9) 2 } . (3.11) 

Z7T 

From the commutation relations it can be seen that the Hamiltonian leads to the equation 
of motion 

d 2 t 9 = v 2 d 2 x 9 

and a similar equation for 0. the velocity v, as well as an additional dimensionless constant 
g depend on the strength of interaction. For non- interacting electrons, v = vf and g = 1. 

It is convenient to expand the original Fermion operator into two part corresponding 
to the motion around the points ±/cp, 

4>{x) w ip R + ipL = e kpX e^ R + e~ kFX e^ m , 

where <&r/l = <P ± 9. These two field commute with one another and satisfy the relations, 

[®r(x), <S>r(x')} = = iir sgn (x - x') . 

The right and left moving electron densities can be reconstructed as 

N B/L = ±d x <S> R/L . 

Then we can rewrite the Hamiltonian (|3.11| ) 



as 



H = nv [N 2 R + Nt + 2\N R N L } (3.12) 



52 



CHAPTER 3. BALLISTIC TRANSPORT 



with 



1 '0 = n I 9 + T I > X = — 2 



2 v <?; i+^ 2 

This Hamiltonian describes interacting system of right and left moving electrons. We 
identify v to vp, the case g < 1 corresponds to A > (repulsion) while p > 1 corresponds 
to attraction. 

It is not straightforward how one can relate the parameters v and g with the original 
lattice model because we have only the effective theory for low-energy excitations. There 
are analytical expressions only for some models, e. g. for the case of long-range interaction 
comparing to the scale hp 1 . 

A very important point is that the parameter g has a physical meaning of dimensionless 
(in units e 2 /h per spin) conductance of an infinite ideal Luttinger liquid. That can been 
traced introducing new variables, <Pr/l = 9<P ± 0, that diagonalize the Hamiltonian ( 3. 10| ) 
as ^ 

H = — (n% + n 2 L ) , n R/L = ±—d x <i) R / L . (3.13) 

The operators Ur/l correspond to the densities of right and left moving interacting elec- 
trons. The Hamiltonian allows chiral general solutions f{x ± vt). Now we can raise the 
chemical potential of right chiral mode ur by an amount fiR Them 57i = —efiRUR, and 
minimizing (|3.13|) we get ur = (ge/2nv)fj,R. Since the additional current to the right is 

Ir = en R v we get 

G = gj. (3.14) 

As is was already mentioned, in a quantum wire it is impossible to couple only to one 
chiral mode. As a result, the d.c. conductance of a finite Luttinger liquid wire coupled 
to noninteracting leads appears the same as for noninteracting case. However, interaction 
effects can be revealed in a.c. experiments, as well as in the presence of some scatterers. 
Because of long-range correlations, the scatterers "undress" the excitations of interacting 
systems. As a result, may interesting and important effects can be observed. In particular, 
the interaction leads to a strong renormalization of density of states near an obstacle. For 
example, if the Luttinger liquid wire has a large barrier with low transmission coefficient 
T one can employ the results for the density of states in a semi-infinite liquid. That results 
in the nonlinear current- volt age curve at low temperatures. 

/ oc T \V\ 2 ^- 1)/9 V -> G(V) = ^oc T \V\ 2{9 - 1)/9 . 

Thus, got the repulsive case the linear conductance is zero. At finite temperature it does 
exist. However it is proportional to T 2 ^ 9 " 1 ^ 9 . It the case of weak scattering the results are 
substantially different. 

There are several experiments where some deviations from the predictions of single- 
electron theory were observed, such as unusual conductance quantization and anomalous 
temperature dependences. Unfortunately, the electron-electron correlations are effectively 
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destroyed by disorder and electron-phonon scattering. Therefore, to observe the interac- 
tion effect one needs extremely pure samples and low temperatures. The results of such 
experiment is demonstrated in Fig. |3.1(J| . 

The concept of Luttinger liquid is specifically important for quantum Hall effect sys- 
tems. We shall see that near the edges of a Hall bar a specific edge states appear which 
can be described by the above mentioned model. This system is much more pure that 
quantum wires, and interaction effects are crucially important. We are going to discuss 
quantum Hall systems later. 
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(a) Wire preparation by cleaved edge overgrowth of GaAs-AlGaAs by MBE. The wire is 
fabricated by cleaving the specimen [see also panel (b)]. Edge states (d) form the quantum 
wire. Panels (e) and (f) show different charge distributions for different top voltages. The 
panel (c) shows a blowup of critical device region. The mean free path is estimated as 10 
/iiii, the length of the channel is about 2 /im. 
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Linear response conductance of a 2 /im log wire in a 25 nm quantum well vs. the top-gate 
voltage (Vt) measured at a temperature 0.3 K. Solid line is the measured conductance. The 
dashed curve is the measured conductance multiplied by an empirical factor 1.15. Inset: 
Linear response conductance of the last plateau for wires of different lengths fabricated 
consecutively along the edge of a single 25 nm cleaved edge overgrowth specimen. The 
numbers denote the wire lengths in microns. 

Figure 3.10: Non-universal Conductance Quantization in Quantum wires [From A. Yacoby, 
et al, Physical Review Letters, 77, 4612 (1996).] 



Chapter 4 

Tunneling and Coulomb blockage 



4.1 Tunneling 



Modern technology allows to fabricate various structures involving tunneling barriers. One 
of the ways is a split-gate structure. Such a system can be considered as a specific example 




Gate 



Figure 4.1: Split-gate structure allowing resonant tunneling. 

of series connection of to obstacles. The complex amplitude of the wave transmitted 
through the whole system is 



D 



1 _ e 2i4> r2r ' l i _ e d ^RJl 2 ' 
where 9 = 20 + arg(r 2 r' 1 ). It is clear that the transmittance 

rp T,T 2 



1 + R 1 R 2 -2y/B^R 2 ~cos9 
is maximal at some specific value of 6 where cos# = 1, the maximal value being 

TiT 2 
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(4.1) 



(4.2) 



(4.3) 
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This expression is specifically simple at Ti,T 2 <C 1, 

4T 1 T 2 



T 

-L 11 



m+T 2 )2- 



(4.4) 



Thus we observe that two low-transparent barriers in series can have a unit transmittance 
if they have the same partial transparencies, T\ = T 2 = T. The reason of this fact in 
quantum interference in the region between the barriers which makes wave functions near 
the barriers very large to overcome low transmittance of each barrier. 

An important point is that the phase 9 gained in the system is a function of the electron 
energy. Thus near a particular value E^ defined by the equality 



cos#(£ (r) ) = -> 9{E i 



(r)v 
k , 



2nk 



one can expand cos 9 as 



2 \dE ) K ' 



Thus at low transmittance we arrive at a very simple formula of a Breit-Wigner type, 

TiT 2 



T 



(T 1 +T 2 ) 2 /4 + (^) 2 (^-^ (r) y 
TjT* 

(r 1 + r 2 ) 2 /4 + (E- J EM) 



2 ' 



(4.5) 



Here we denote 9' = (d9/dE) E _ E ( r ) and introduce = Ti/\9'\. 

The physical meaning of the quantities Tj is transparent. Let us assume that all the 
phase shift is due to ballistic motion of an electron between the barriers. Then, 



9 = 2ka = 2ah'W2mE -> 9' 



2m 
~E 



2a 
hv 



frv a Ti , 



where v is the electron velocity. As a result, the quantity Tj can be rewritten as T 
where v a = v/2d is the frequency of oscillations inside the inter-barrier region, the so-called 
attempt frequency. Thus it is clear that Tj are the escape rates through z-th barrier. 

To specify the transition amplitudes let us consider a ID model for a particle in a well 
between two barriers. We are interested in the scattering problem shown in Fig. |4.2| . To 
find a transmission probability one has to match the wave functions and their gradients at 
the interfaces 1-4 between the regions A-C. They have the following form 
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Figure 4.2: On the resonant tunneling in a double-barrier structure. 



Here 

k = TC x \llmE , Ki — \J /«Qj — k 2 , k i = FT 1 \f2mUi . 

The transmission amplitude is given by the quantity t while the reflection amplitude - by 
the quantity r. In fact we have 8 equations for 8 unknowns (r, t, a«, bi, q), so after some 
tedious algebra we can find everything. For a single barrier one would get 

Ak 2 K 2 k 2 K 2 

T(E) = 4/c K ~ K 1 c -"2Kd 

y ! 4sm\i 2 {Kd)+Ak 2 K 2 ~ 4 

Here d is the barrier's thickness. So the transparency exponentially decays with increase of 
the product nd. The calculations for a double-barrier structure is tedious, so we consider 
a simplified model of the potential 

U(x) = U d[5(x) + 5(x - a)} . 

In this case we have 3 regions, 

e ikx _|_ r( ,-ikx X < . 

A sin kx + B cos kx < x < a , (4.7) 

te ik(x-a) x> a 

The matching conditions for the derivatives at the 5-functional barrier has the form 

tf(xo + 0) - rj}\x Q - 0) = K 2 d^j(x ) . (4.8) 

Here k 2 = 2mU /H 2 . One can prove it by integration of the Schrodinger equation 

(h 2 /2m)V 2 ^ + U d5(x - x )ip = 

around the point x . Thus we get the following matching conditions 

B = 1 + r, 

kA — ik(l — r) = /t 2 a(l + r), 

A sin ka + B cos ka = t , 

ikt — k(A cos ka — B sin ka) = tn 2 a. 
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First one can easily see that there is a solution with zero reflectance, r = 0. Substituting 
r = we get the following requirement for the set of equation to be consistent 

k = k , tan k Q a = — . (4.9) 

We immediately observe that at that k \t\ = 1 (total transmission). At strong enough 
barrier, nd 3> 1, this condition means 

k a = 7r(2s + 1) , s = 0, ±1, .. 

Physically, that means that an electron gains the phase 2ns during its round trip (cf. with 
optical interferometer). Thus two barriers in series can have perfect transparency even if 
the transparency of a single barrier is exponentially small. The physical reason is quantum 
interference. 



h 2 kl 

2m 



The condition (|4.9|) defines the energy 

En 



where the transparency is maximal. Near the peak one can expand all the quantities in 
powers of 

E — Eq E — Eq 



k — k « : « fu - 



(dE/dk) ko 2E 
The result for a general case can be expressed in the Breit-Wigner form 

r L r R 



T{E) 



(E-E Q y + \{T L + v R y 



Here Tur^/H are the escape rates for the electron inside the well to the left (right) lead. 
They are given by the attempt frequency Vn/2a = hkn/2ma times the transparency of a 
given barrier. 

Of course, if voltage across the system is zero the total number of electrons passing along 
opposite directions is the same, and the current is absent. However, in a biased system 
we obtain the situation shown in Fig. [4.3| . Negative differential conductance, dJ/dV < 0, 
allows one to make a generator. One can also control the system moving the level Eq with 
respect to the Fermi level by the gate voltage. In this way, one can make a transistor. 

Commercial resonant tunneling transistors are more complicated than the principle 
scheme described above. A schematic diagram of a real device is shown in Fig. |4.4| . In this 
device resonant tunneling takes place through localized states in the barrier. There exist 
also transistors with two quantum wells where electrons pass through the resonant levels 
in two quantum wells from the emitter to collector if the levels are aligned. The condition 
of alignment is controlled by the collector-base voltage, while the number of electrons from 
emitter is controlled by the base-emitter voltage. 
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■ >>■ eV 

^X7 



Bias voltage 

Figure 4.3: Negative differential conductance in double-barrier resonant-tunneling struc- 
ture. 



4.2 Coulomb blockade 

Now let us discuss a specific role of Coulomb interaction in a mesoscopic system. Consider 
a system with a dot created by a split-gate system (see above). 

If one transfers the charge Q from the source to the grain the change in the energy of 
the system is 

AE = QV G +^. 

Here the first item is the work by the source of the gate voltage while the second one is the 
energy of Coulomb repulsion at the grain. We describe it by the effective capacitance C to 
take into account polarization of the electrodes. The graph of this function is the parabola 
with the minimum at 

Q = Q = -CV G , 



So it can be tuned by the gate voltage V G . 
transferred by the electrons with the charge 
number n of electrons at the grain is 



Now let us remember that the charge is 
— e. Then, the energy as a function of the 



AE(n) = -neV G + — . 
Now let us estimate the difference 

e 2 

AE(n + 1) - AE(n) = -eV G + n— . 

C 

We observe that at certain values of V G , 

V Gn = n^, (4.10) 

the difference vanishes. It means that only at that values of the gate voltage resonant 
transfer is possible. Otherwise one has to pay for the transfer that means that only inelastic 
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Figure 4.4: Schematic diagram of a Si MOSFET with a split gate (a), which creates a 
potential barrier in the inversion layer (b). In the right panel oscillations in the conductance 
as a function of gate voltage at 0.5 K are shown. They are attributed to resonant tunneling 
through localized states in the barrier. A second trace is shown for a magnetic field of 6 
T. From T. E. Kopley et al, Phys. Rev. Lett. 61, 1654 (1988). 



processes can contribute. As a result, at 



the linear conductance is exponentially small if the condition (|4.10 ) is met. This phe- 
nomenon is called the Coulomb blockade of conductance. 

As a result of the Coulomb blockade, electron tunnel one-by-one, and the conductance 
vs. gate voltage dependence is a set of sharp peaks. That fact allows one to create a 
so-called single- electron transistor (SET) which is now the most sensitive electrometer. 
Such a device (as was recently demonstrated) can work at room temperature provided the 
capacitance (size!) is small enough. 

Coulomb blockade as a physical phenomenon has been predicted by Kulik and Shekhter [ED] 
There are very good reviews [[H, EH, 113] about single-change effects which cover both prin- 
cipal and applied aspects. Below we shall review the simplest variant of the theory, so 
called "orthodox model" . 



A simple theory of single charge tunneling 

For simplicity, let us ignore discrete character of energy spectrum of the grain and assume 
that its state is fully characterized by the number n of excess electrons with respect to 
an electrically neutral situation. To calculate the energy of the systems let us employ 
the equivalent circuit shown in Fig. [4.5| . The left (emitter) and right (collector) tunnel 



junctions are modeled by partial resistances and capacitances. 
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Figure 4.5: Equivalent circuit for a single-electron transistor. The gate voltage, V g , is 
coupled to the grain via the gate capacitance, C g . The voltages V e and V c of emitter and 
collector are counted from the ground. 



The charge conservation requires that 

— Ue = Qe + Qc + Qg 

= C e (V e -U)+C c (V c -U)+C g {V g -U), (4.11) 
where U is the potential of the grain. The effective charge of the grain is hence 

Q = CU = ne+ CM > C = J2 c i- 

i=e,c,g i 

This charge consists of 4 contributions, the charge of excess electrons and the charges 
induced by the electrodes. Thus, the electrostatic energy of the grain is 



p Q 2 inef 
n 2C 2C 



+ ^E^+^(e^) • ( 4 - 12 ) 



The last item is not important because it is n-independent. In the stationary case, the 
currents through both junctions are the same. Here we shall concentrate on this case. In 
the non-stationary situation, an electric charge can be accumulated at the grain, and the 
currents are different, see Appendix 0. 

To organize a transport of one electron one has to transfer it first from emitter to grain 
and then from grain to collector. The energy cost for the first transition, 

E n+ i-E n = ^^ + ^ C M ( 4 - 13 ) 

i 

must be less than the voltage drop eV e . In this way we come to the criterion 

E n - E n+1 + eV e > . (4.14) 

In a similar way, to organize the transport from grain to collector we need 

E n+1 -E n -eV c >0. (4.15) 
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The inequalities (|4.15|) and (|4.14|) provide the relations between V e , V c and V g to make the 
current possible. For simplicity let us consider a symmetric system, where 

G e = G c = G, C e = C c ^ C/2 {C g < C), V e = -V c = V b /2 

where V b is bias voltage. Then we get the criterion, 

V b > {2n + l)\e\/C-2{C g /C)V g . 

We observe that there is a threshold voltage which is necessary to exceed to organize 
transport. This is a manifestation of Coulomb blockade. It is important that the threshold 
linearly depends on the gate voltage which makes it possible to create a transistor. Of 
course, the above considerations are applicable at zero temperature. 
The current through the emitter-grain transition we get 



J2pn[r e - >g -r g ^ e } . (4.16) 



Here p n is the stationary probability to find n excess electrons at the grain. It can be 
determined from the balance equation, 

Vn-iK-i + Pn + iK +1 - (rr 1 + r n n+1 ) Pn = o. (4.17) 

Here 

r n n _, = r e _ s (n-l) + r c _ s (n-l); (4.18) 

= T g ^ e (n + 1) + T g ^ c (n + 1). (4.19) 

The proper tunneling rates can be calculated from the golden rule expressions using tun- 
neling transmittance as perturbations. To do that, let us write down the Hamiltonian 

as 



e q C qcr C qo- ) 



H 9 = J2' 

qcr 

U ch = (n-Q )/2C, n = Y, C i« c w- N " 



q<r 



Here 7i bath is the Hamiltonian for the thermal bath. We assume that emitter and collector 
electrodes can have different chemical potentials. iV + is the number of positively charged 
ions in the grain. To describe tunneling we introduce the tunneling Hamiltonian between, 
say, emitter and grain as 

T^-e^g = ^ ^kqCjtcrCqfj + h.C. 
k,q,(T 



4.2. COULOMB BLOCKADE 



63 



Applying the golden rule we obtain 

, G 



/oo poo 
de k / de q / e (e fe )[l - f g (e q )] 5(E n+1 - E n - eV e ) . 
■oo J — oo 



Here we have introduced the tunneling conductance of e — g junction as 

G e = (47ce 2 /h) g e (e F )g g (e F ) V e V a (|T kq | 2 > . 

along the Landauer formula, V ej9 being the volumes of the lead and grain, respectively. In 
this way one arrives at the expressions 

2a 



Here 



while 



r e _> g (n,V e ) 



" g ^c(n, V c ) 



(~n, -V e ) 



T c ^ g (-n, -V c 



+,e) j 



(4.20) 
(4.21) 



1 + exp(-e/fcT) 
1 



A±,»(n) = E n - E n±1 ± eV, 



C 



e0(e) at T -> 



=F en =f e C,^ 



is the energy cost of transition. The temperature-dependent factor arise from the Fermi 
occupation factor for the initial and final states, physically they describe thermal activation 
over Coulomb barrier. The results of calculation of current-voltage curves for a symmetric 
transistor structure are shown in Fig. |4. At low temperatures and low bias voltages, 




Figure 4.6: The current of a symmetric transistor as a function of gate and bias voltage at 
T = (from the book [0]). 



VC /e < 1, only two charge states play a role. At larger bias voltage, more charge states are 
involved. To illustrate this fact, a similar plot is made for symmetrically biased transistor, 
V e = -V g = V/2, for different values of Q , Fig. [47. 
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Figure 4.7: The current of asymmetric transistor, G e = 10G C , as a function of bias voltage 
at T = and different Qo-e = 0, 0.25 and 1 (from the book [0]). At Qo = the Coulomb 
blockade is pronounced, while at Qo/e = 0.5 the current-voltage curve is linear at small 
bias voltage. The curves of such type are called the Coulomb staircase. 

Cotunneling processes 

As we have seen, at low temperature the sequential tunneling can be exponentially sup- 
pressed by Coulomb blockade. In this case, a higher-order tunneling process transferring 
electron charge coherently through two junctions can take place. For such processes the 
excess electron charge at the grain exists only virtually. 

A standard next-order perturbation theory yields the rate 



Two features are important. 

• There are 2 channel which add coherently: (i) e — > g, g — > c with the energy cost 
A_ >e (n + 1), and (ii) g — ► c, e — > g with the energy cost A +jC (n — 1). 

• The leads have macroscopic number of electrons. Therefore, with overwhelming 
probability the outgoing electron will come from a different state than the one which 
the incoming electron occupies. Hence, after the process an electron-hole excitation 
is left in the grain. 

Transitions involving different excitations are added incoherently, the result being 




2 



r. 



cot 




Je Jg Jg Jc 



/ de k / de q / de q , \ de k , f(e k )[l - f(e q ))f(e q r)[l - f(e k ,)] 



4.2. COULOMB BLOCKADE 



65 



At T = the integrals can be done explicitly, and one obtains 



r 



hG e G c r i i 



V 3 for eV < Aj . 



cot 



127re |_A_ e (n + l) A +iC (ra - 1) 



As a result, the current appears proportional to V 3 that was observed experimentally. The 
situation is not that simple for the degenerate case when Aj = 0. In that case the integrals 
are divergent and the divergence must be removed by a finite life time of a state. A detailed 
treatment of that case is presented in the book [0] . 

There is also a process when an electron tunnels through the system leaving no excita- 
tions in the grain. The probability of such elastic cotunneling has a small factor (g^Vg) -1 . 
However, it leads to the current, proportional to V, thus it can be important at very low 
bias voltage. 

Concluding remarks 

There are many experiments where Coulomb-blockaded devices are investigated. Probably 
most interesting are the devices where tunneling takes place through a small quantum 
dot with discrete spectrum. An example of such device is shown in Fig. The linear 
conductance of such a structure as a function of the gate electrode C is shown in Fig. fO| 
An important point is that at present time the devices can be fabricated so small that the 
criterion kT < e 2 /C can be satisfied at room temperatures. Now several room temperature 
operating Coulomb blockade devices were reported. Among them are devices consisting 
of large molecules with the probes attached in different ways. This is probably a starting 
point for new extremely important field - molecular electronics. Such devices are extremely 
promising both because they are able to operate at room temperatures and because they 
will allow high integration. This is one of important trends. Another one concerns with 
single-electron devices which include superconducting parts. There is a lot of interesting 
physics regarding transport in such systems. 
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Figure 4.8: (a) A typical structure of a quantum dot. The depleted (shaded) areas are 
controlled by electrodes 1-4, C, and F. Electrode C also controls the electrostatic potential 
in the dot. (b) a model of a quantum dot. From [□]. 




0.50 



Figure 4.9: Conductance of a quantum dot vs. the voltage of gate electrode C. From L. P. 
Kouwenhoven et al, Z. Phys. B 85, 367 (1991). 



Chapter 5 



Conductance Fluctuations. 
Mesoscopics 



5.1 General considerations 



We have mentioned that fluctuations between the values of conductance of different samples 
become very large in low-dimensional systems. However, the samples still contain many 
atoms. According to the estimates, one can expect a very specific behavior in the sample 
with the size ~ 1(T 4 cm at T = 0.1 K. So, 

• where is the difference between macroscopic and microscopic systems? 

• What is going on in between, in such a so-called mesoscopic region? 

It is very important that one can observe mesoscopic behavior even in a given sample 
monitoring its properties as a function of an external parameter, say magnetic field. The 
reason is that these effects are manifestations of quantum interference. External magnetic 
field effects the phase gains along the trajectories which behave as almost random quanti- 
ties. So (with some restrictions) one can consider a sample in a given field as a separate 
representation of the ensemble of samples. Experimentally, mesoscopic fluctuations man- 
ifest themselves as fluctuations in the resistance as a function of magnetic field. They 
look random. However, they are "fingerprints" of a sample, they are reproduced at any 
experimental run. 

Now we discuss mesoscopic fluctuations in more detail. 

5.2 Universal conductance fluctuations 

As we have already mentioned, in 2D case the Drude conductance for a single spin direction 
(and a single valley) can be written as 



G = 



~h2L 



N , N = 



k F W 



7T 
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Here L is the length of the sample, W is its width, while £ is the elastic mean free path. 
The number N is just the number of subbands that are occupied at the Fermi energy in a 
conductor with width W. 

Let us have in mind the Landauer picture of two reservoirs with some disordered region 
in between.^] For simplicity let us assume that the latter region is connected to reservoirs 
by ideal quantum wires. 

If L > £ then the modes have the same average transmission probability, ir£/2L, 
that we can establish comparing Drude and Landauer formulas. We are interested in the 
fluctuations around this average. For that let us come back to the multichannel Landauer 
formula (for one spin) 

e 2 N 

G = -^ Yl i^i 2 - 

a,/3=l 

For the ensemble average transmission coefficient we get 

(\t aP \ 2 ) = 7c£/2NL. 



We are interested in the quantity 

Var(C7) = <(C7-(G» 2 ) . 

An important point that while doing such an average one should keep track of the correla- 
tion in transmission probabilities |t a ,g| 2 for different pairs of incident and outgoing channels. 
The reason is that the transmission takes place via many scattering events, and they may 
involve the same scatterers. On the other hand, the reflections are dominated only by 
few scattering events, so they can be considered as uncorrelated. These considerations 
definitely need a rigorous derivation which has been given first by Altshuler [EEl]. Taking 
this concept for granted, we get 

N N 

M 2 = n- £ M 2 (5.i) 

a, /3=1 a,/3=l 

so the variance of the conductance is 



Var(G)= (£) Var \JT, |r Q/3 | 2 j 



2\ 2 

Ti 



iV 2 Var (\r a/3 \ 2 ) , (5.2) 



assuming reflections to be uncorrelated. There are different Feynman paths for reflection, 
so we can write the variance 

Var (|r Q/3 | 2 ) = <|r Q/3 | 4 ) - <|r a/3 | 2 ) 2 
1 Here we follow the arguments by Lee [27]. 
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using summation over the paths, 

M 

<M 4 ) = £ {A*(i)A{j)A*{k)A{l)) 

i,j,k,l=l 
M 

= E [<i^(i)i 2 ><i^(*)r)M« + (i^wr)(i^(j)i 2 )^] 

i ,j,k,l=l 

= 2{\r aP \ 2 ) 2 . (5.3) 
Here we neglected the terms smaller by a factor 1/M <C 1 . Therefore 

Var (|r Q/3 | 2 ) = (|r a/3 | 2 ) 2 . 



Since the average 

<M 2 > ' 



AT 



1 - O 



L 



that follows from the averaging of ( |5.1|) , we come to the conclusion that Var (G) = (e 2 /h) 2 
independently of £ and L in the diffusive limit f C L. This is why the called the universal 
conductance fluctuations (UCF). 
The general formula reads as 

5G = VVMG) = ^LC^- . 
v v ; 2V3 h 

Here gi are degeneracy factors, C depends on the sample geometry (actually, of its effective 
dimensionality), while j3 = 1 in zero magnetic field and 2 when the magnetic field breaks 
time- reversal symmetry (see later). 

5.3 Temperature effects 

The thermal effects are a bit complicated because of finite phase coherence length = 
y^Dr 9 and of thermal averaging characterized by the thermal length = \JhD/kT. Both 
effects tend to restore self-averaging. We discuss the situate in a semi-qualitative way for 
the geometry of a narrow channel in a 2DEG, namely at 

W < < L . 

At L 9 Lt one can neglect the effects of thermal average and think only about coherence. 
Then one can divide the sample into un correlated segments with the length L v , each of 
which has a fluctuation of order e 2 /h. All the segments are in series, so their resistances 
should be added according to the Ohm's law. Denoting the resistance of each segment by 
Rtp we get 

Var {R v ) w (i^) 4 Var (!/#„) = (R^) A (e 2 /h) 2 . 
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The average resistance is (R) = (L/L ip )(R lfi ). Thus the variance of the resistance is 

Var (R) = (L/^)Var(i^) = {L/ L v )(R v )\e 2 /h) 2 . 
Coming back to conductance, we get 

Var(G) = (R) 4 Var(R) « (L v /L) 3 (e 2 /h) 2 . 

Finally 

5G = constant x j . (5.4) 

Now let us discuss the role of thermal average. Consider two Feynman paths with the 
energies different by 5E. The can be considered as uncorrelated after a time t\ if the 
acquired phase difference Sip = (5E)ti/h is of the order unity. In this time the electron 
diffuses a distance L\ = \/Dt\ rs yhD/ (8E). Now we can define the correlation energy at 
which the phase difference following diffusion at the distance Li is unity, 

E^^HD/L 2 . (5.5) 

Equating E c = kT we obtain the expression for L T . If ^> L T , then the total energy 
interval kT near the Fermi level is divided into subintervals of the width E C (L V ) = h/r^ in 
which the phase coherence is maintained. These intervals can be assume as uncorrelated, 
f] their number being N ~ kT/E c (L v ). As a result, the fluctuations are suppressed by 
the factor N~ 1 / 2 = L T jL^. There exists a very useful interpolation formula (see, e. g. 
Ref. [ED]) 

-1/2 



5G 



g s g v \/l2 e 2 ( " J 



2^ h V L 



9 fL v 



2vr \L t 



(5.6) 



5.4 Magnetoconductance fluctuations 

In real experiments, people change either Fermi energy (by gates) or magnetic field H 
rather than impurity configuration. It is conventionally assumed that sufficiently large 
changes in the Fermi energy or in magnetic field is equivalent to a complete change of 
impurity configuration (the so-called ergodic hypothesis"). The reason is that there exist 
characteristic values of the correlation energy E c or magnetic field, H c . That provides a 
possibility to study ensemble averages without changing the sample. The magnetoconduc- 
tance correlation function is defined as 

F(AH) = ({5G{H) - (G(H))] [5G{H + AH) - (G(H + AH))]) . 

Thus the previous correlation function is just F(0) while the correlation field H c is defined 
as F(H C ) = F(0)/2. 

To understand what is going on let us take into account that the correlation func- 
tion contains the product of 4 Feynman path amplitudes, A(i, H), A*(j, H), A(k, H + 
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H + AH 






(e/h) (<£ + A<£) 



Figure 5.1: Illustration of the influence of magnetic field on the diffuson (left) and cooperon 
(right) contributions. 



AH),A*(l,H + AH) along various paths between the points r and r', see Fig. 5.1 . The 
term for which i = I, j = k, i. e. 

A(i, H)A*(j, H) x A*(i, H + AH)A(j, H + AH) 

is called the diffuson. The physical meaning of this term is comes from the understanding 
that A(i, H), A*(i, H) is just the classical probability for a diffusion from the point r to 
the point r'. The phase of this term in a magnetic field is 



A • d\ + - 

el c 



A + AA) ■ d\ 



h 



-A$. 



It contains only the increment of the flux, A$ = A ■ AH and does not contain the flux 
A ■ H itself. 

For the cooperon which is the product 

A. (i, H)A*_ (j, H)A+(j, H + AH)A* + (H + AH) 

where the — sign refers to the trajectory from r' to r (i. e. time-reverted) while + sign 
refers to the trajectory from r to r', the phase is 



A ■ d\ + - 

c I c 



'A + AA) ■ d\ 



h 



(2$ + A$) 



In contrast to the diffuson, the cooperon is sensitive to the total flux $ through the loop 
and therefore can be suppressed by a weak magnetic field. 

We conclude that at B > B c only the diffuson contributes to the magnetoconductance 
fluctuations. The combined effect of magnetic field and inelastic scattering on the diffuson 
can be allowed for by the replacement r v — > r cfi with 



'off V 1 "Aff/2 • 

Here t ah/2 corresponds to the magnetic relaxation time tn obtained for weak localization 
with substitution H — > AH/2. 

2 This is actually true only at W L v . 
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5.5 / — V-curve of a metallic wire: Fluctuations 

Let us consider an example - current-voltage curve of a small (but macroscopic) metallic 
wire with 

£ < b < L < . 

Let us estimate the phase difference between the electrons which have the difference in 
energies 5e and had a total path S within the sample: 

S\p(e + 6e) - p{e)\ S 5e 
h hv 

The electronic states can interfere if A(p c < 1. On the other hand, for a diffusive motion 
S ~ vt ~ vL 2 /D. Thus 

hv hD 
S L 2 

As a result, I — V curve must fluctuate at the voltage scale 

V c ~ hD/(eL 2 ) . 

What is the scale of the fluctuations in current? To estimate it let us take into account 
that only the electrons within the layer eV near the Fermi level are important (at T = 0). 
However, the fluctuations within the strips having the width V c are statistically indepen- 
dent, the number of such intervals being V/V c (at V ^> V c ). Each interval contributes to 
the current as 

VA . / V r \ V 2 V c e 2 

b (Acr) ~ 



RJ \L/ab 2 J L h ' 

Here we have assumed Id formula Aab 2 ~ (e 2 /h)L. Making use of statistical independence 
of fluctuations, we obtain 

a 

Here we used the similarity to the diffusive motion: 

x ~ \Tb~t ~ Vvti~ Vse 



with the mapping 
It is interesting that 



c ■ 



h e 2 VW c jv 

V^h V n ~ hVVr. 



' c '° ' c "-V^c 

It increases with V. This one can face negative differential conductance and corresponding 
instabilities at 

v ( h 

> 



V c V e 2 R c 
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As we have seen, the interference interval is crucially important. At finite temperature, 
the energy spread T plays the same role as the interval eV above. Consequently, we have 
T I {be) independent intervals, the summary effect being proportional to \/T5e. Thus the 
average manifestation of the interference goes as 

VT5e [Se 
T V T ' 

Requesting 5e > T we get 

L < Lt ~ y/hD/T . 

Thus in general (with some restrictions) one has to require L < L v ,Lt to observe meso- 
scopic phenomena. 



5.6 Random matrix theory 

Universal conductance fluctuations (UCF) can be described by a very different approach 
based on on statistical theory of energy levels in complex systems, developed originally in 
atomic physics (see [DZZI] for a review). The approach was initially suggested by Altshuler 
and Shklovskii ['29] and Imry [3(1]. The starting point is the Thouless formula for the 
dimensionless conductance Q = G/Gq, Go = e 2 //i, 

G = E C /A, (5.7) 

where E c is a typical shift of the discrete levels in a finite sample induced by changing the 
boundary conditions, say, from periodic to antiperiodic, while A = [g(eF)L d ]~ l is a typical 
interlevel spacing at the Fermi energy. 

The physical meaning of Eq. ( |5.7|) can been understood having in mind that E c is just 
the quantum mechanical uncertainty associated with the time required to move through 
the sample. Thus, for the diffusive system, 

E c = hD/L 2 , (5.8) 

where D is the diffusion constant. By the contrast, it the states are localized, E c decays 
exponentially with the sample length L, E c ~ exp(— L/£) where £ is the localization length. 
If we formally introduce the scale- dependent diffusion constant, D(L), then we come from 
Eq. (|5.7|) to a generalized Einstein relation, 

Q{L) ~ hg{e F )D{L)L d - 2 

which is just the Drude formula for D = (l/d)vpL 

According to Eq. (|5.7| ), the reduced conductance Q equals to the number of levels, N c , 
within the band of width E c at the Fermi level. In metallic regime, N c ^> 1. One could 
expect that that the variance of the quantity N c defined as 



y ai N c =(N 2 )-(N c ) 2 , 
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is of the order N c . However from the microscopic theory it follows that 

VarJV c «l. (5.9) 

From this point of view, the fluctuations are anomalously small, though they are much 
larger than the thermodynamical ones. 

The basic reason for that is that interlevel spacings are correlated. To quantitative 
treatment of the correlation the Wigner- Dyson approach is usually used. The approach 
assumes that the energy eigenvalues of complex systems may me approximated by those 
of random matrices 7i for their Hamiltonians. 

For pedagogical reasons, let us first discuss the case 2x2 symmetric random matrices, 

n _ ( hu h\2 

~ V ^12 h 22 

with 3 independent random parameters hu, h 12 , h 22 . The eigenvalues are 

E± = {h ± ±s)/2; 
h± = h lx ± h 22 , 



s = Jh 2 _+Ah 2 



12 



Let us find the distribution of the eigenvalues, p{E + , E-), assuming known the distributions 
of random matrix elements, V{h\\, h\ 2 , h 22 ). We have 

p(E + ,E^) = J dh 11 dh 2 2dh 11 V{h 11 ,h 12 ,h 2 2) 

xS(E + - h + /2 - s/2) 5 {E_ - h + /2 + s/2) . 

Since 

dhu dh 22 = 2dh + dh_ and dh + ds = dE + dE_ 

we get 

s f s dh_ n(E + -E_) 
oc . 



P (E + ,E-)«Jjh. -g- J_ i7 = 



Thus we conclude that the probability to find two eigenvalues vanishes as \E 2 — Ei\ for 
close levels. This is the so-called quantum level repulsion. 

The original Wigner-Dyson (WD) prescription which allows to prove several important 
theorems for the so-called Gaussian orthogonal ensemble is 

p(H) oc exp (-^Tr(H 2 ) 

where N is the matrix range while A is a numerical parameter. Using this assumption we 
obtain 

(J? tc(E + - EJ) ( 2{E\ + El 
P{E+,E_) = -± '- exp ' 



\ 2 
A o 
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The probability for the interlevel spacing s is thus 

1 f°° ( s 2 

P(s) = g J dh + p(E + ,E^) oc s exp f 

This is the famous WD surmise for the spacing distribution. It is exact for 2x2 matrices, 
but also tends to be true for any N. Indeed, for 2 very close levels the problem is effectively 
reduced to a 2 x 2 matrix, and we arrive at the expression 



p(s) = J dh„ J dh 12 S[s- \Jh 2 _ +lh\^j 



which is proportional to s at small s. 

Now let us turn to the case when the time-reversal symmetry is broken (say, by an 
external magnetic field). Then the corresponding Hamiltonian is complex (and Hermitian) 
rather than real symmetric, and the level spacing is 



h 2 _ + (2^h 12 ) 2 + (2^h 12 y. 

Now three parameters must vanish simultaneously for an accidental degeneracy s = 0. As 
a result, at small s 

p(s) oc s 2 . 

The general expression for the eigenvalue distribution can be constructed as 



p(Ei <E 2 ...<E n )(x Y[(Ei - Ejf exp ( V 

i>3 V A ° - 



E 



Here /3 — 1, 2 and 4 for the so-called orthogonal, unitary and symplectic ensembles. 

To prove the adequacy of the description one should calculate the conductance variance 
and compare it with microscopic derivation. To calculate the variance the Dyson-Mehta 
theorem can be used. It states that the variance of any random variable of the form 

A = Y,<E l ) 

i 

can be expressed as 



Var A = -— / tdt 

n 2 l3 Jo 



J dEe iEt a(E) 



For our purpose of evaluating N c we assume 

a{E) = 1 - Q(\E\ - E c /2) . (5.10) 



3 A symplectic matrix is a unitary matrix with real quaternion elements. 

4 This expression is called the linear statistics because it does not intermix different eigenvalues. 
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Applying the Dyson-Mehta theorem we find that the expression for VariV c diverges at 
small t unless the integral is cut off at t ~ A -1 . As a result, 

This formula shows strong depression of fluctuations. However it differs from the micro- 
scopic derivation by a large logarithm. The discrepancy was analyzed by Altshuler and 
Shklovskii who pointed out that the above expression is actually valid for a closed sys- 
tem. In an open system the levels are broadened by the amount of the order E c itself. To 
incorporate this broadening one can replace Eq. ( |5.10|) by a smooth-edge expression 



a(E) = (1 + E 2 /E, 



to obtain Var N c = 1/4/3. Unfortunately, the numerical factor here depends on the partic- 
ular shape of a(E) and not universal. A generalization which yields the exact numerical 
factors exists only in ID case (see [Id] for a review). This is the case of 2-terminal resistance 
between the leads 1 and 2 which is given by the Landauer formula 

£ = 2TrS 21 4i- 

The matrix S^iS"^ is Hermitian and positive, it has N transverse eigenvalues corresponding 
to N transverse modes with < Tj < 1. The result for the average distribution function 
for the parameters Aj = (1 — Tj)/Tj has the form 

P (A) = 2M ' L , 
VA(1 + A)' 



where it is assumed that N£ 3> L ^> I. 



Chapter 6 

Quantum Hall Effect and Adiabatic 
Transport 

6.1 Ordinary Hall effect 

In a magnetic field electrons' trajectories are curved because of Lorentz force. As a result, 

j = a (E + [j x H]/nec) , a = ne 2 T/m . 
One can solve this vector equation to obtain the resistivity tensor 

, Po = l/°o • 



P = 



Po H/enc 
-H/enc po 



The inversion of this tensor gives the conductivity tensor with the components (in 2d case) 

_ nee 1 
1 + (u c t) 2 H uj c t 

There is a striking similarity between the quantization of the conductance of a ballistic 
channel in units of e 2 /h and of the Hall conductance. 

6.2 Integer Quantum Hall effect - General Picture 

In the quantum case one faces the Landau levels. We have seen that the number of states 
per unit area for a filled Landau level is 

uh = l/2na 2 H = eH/ch . 

Usually the filling factor 

v = n/nn 

for a fractionally filled level is introduced. If one substitutes to the expression for the Hall 
component of the conductivity tensor and assumes uj c t — > oo he obtains a xy = ue 2 /h. 
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This result seems very fundamental. Indeed, according to the electrodynamics, one can 
eliminate the electric field on a spatially-homogeneous system by assuming a drift velocity 

v = c[E x H]/H 2 . 

Thus, the result seems fundamental independently whether classical of quantum mechanics 
is discussed. 

Experimentally, the following picture was observed. It is clear that only violation of 
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Figure 6.1: Schematic dependence of Hall resistance on filing factor. 

the translational invariance can lead to such a picture. Thus one has to consider either 
impurities, or edges. 

The generally accepted picture is as follows. Impurities violate the translational invari- 
ance. As a result, p y is not a good quantum number any more, and Landau levels smear 
into sub-bands. The crucial point that the most of the states are localized and cannot carry 
the current. 

To make the analysis as simple as possible let us discuss a 2d electrons in crossed electric 
and magnetic fields (E || x, H || z) and a single impurity at the point r = {xo,yo}. As 
we have seen (see Shubnikov-de Haas effect), the weak electric field leads to the energy 
shift p y v where v = cE/H is the drift velocity in ^-direction, as well as to the shift in 
the center-of-motion co-ordinate is shifted in x-direction by v/uj c . Using the corresponding 
states as a basis, we can now expand the exact wave function as 



E 

np y 



C-npy^npy ( r ) 



We get 



E' 

npy 



■np v 



n + v 



^Pnpy Cnp y [E n p y + V] Ipnpy E C n p y 1p np , 



npy 



npy 



6.2. INTEGER QUANTUM HALL EFFECT - GENERAL PICTURE 



79 



Because V oc 8(r — r ) one can write 



A 



E-E. 



np y 



Now we recall that J2 nPy c nPy ip nPy (r ) = ^(r ). Substituting the previous expression into 
this equation, we get the exact condition for eigen energy 



A ^ 



\4>n,py(ro)\' 

E-E, 



n,p y 



The right hand side of this equation as a function of the energy is shown below One can 
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Figure 6.2: Formation of localized states in 2DEG in magnetic field. 

find from this equation one completely localized state for each Landau level, its energy 
shift being proportional to A. The lowest such state can be represented as (for r = 0) 



V>ioc ~ exp 



ixy 
2^ 



x 2 + y 2 
4a 2 H 



The other levels are almost unperturbed and extended. 

Now let us take into account only the lowest Landau level which we assume to be 
completely filled. We have one localized state and N—l extended ones where iV = A/2na 2 H , 
which we can label by the discrete quantum number k as 



k = 



Py^y 
2ith 



Each mode behaves just as the transverse mode in a quantum channel, and the current is 
given as 

I = ~J L ^2( E n,k+i - E nk ) — —- ^(-E„,fc max - E n>kmin ) . 

nk nk 

It is not trivial to prove this equation. It was done by R. Prange using gauge considerations. 
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Proof : 

The main procedure is a s follows. We have specified periodic boundary conditions 
along y-axis. Consider the system as a cylinder and introduce an auxiliary constant 
vector-potential along y axis as 2ira/L y . The vector-potential of such a type can be 
eliminated by the gauge transform 

ip — > exp(2nay / L y )ifj . 

This function can satisfy periodic boundary conditions only if a is integer. Thus the 
extended states which extend from to L y must depend on a. 

The current operator can be written as I — (—e/h)ctH/da, while the average current 
is 

nk 

According to the construction of the quantum number k the introduction of the 
vector-potential leads to the replacement k — > k + a. Thus 

E n ,k\a=l — E Hi k + i\ a= Q . 

Replacing the derivative by the average value we get the result given above. 
Thus we come to the following picture, as the Fermi level passes the regions of the 
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Figure 6.3: Density of states in 2DEG in magnetic field. 

extended states the steps in Hall resistance and peaks at the longitudinal resistance occur. 
As we have shown, the current is independent of the density of states, only the number of 
occupied extended states is important. 

Now we have to remember that the state with A; max (^min) correspond to the upper 
(lower) edge of the sample if we map the quantum number k to the centers of gravity of 
the states. Thus we come in a natural way to edge states. 
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6.3 Edge Channels and Adiabatic Transport 

The quantization of the conductance of the ballistic channel arises from the finite number of 
propagating modes each of which can carry only a very specific current. Thus it is tempting 
to consider the modes of an ideal finite system embedded into an external magnetic field. In 
this simplified picture we can obtain some understanding concerning the nature of localized 
and extended states. 

Let us start from a ideal electron system confined by the potential V(x) in the presence 
of the magnetic field H || z. For a single spin component we have the Hamiltonian 

It is natural to look for a solution in the form (7i commutes with p y \) 

K k) =^nk(x)e iky 

where hk is the eigenvalue of p y . The average velocity for such a state is 



v n (k) = (n,k 



p y + (eH/c)x 



n,k) = (n,k 



on 



d Py 



ldE n {k) 



m 

It is easy to calculate such a velocity for a parabolic confinement, 

V(x) = -rauo^x 2 . 

The result is 

, . hk hk 1 

V{ } = M = m 1 + (ujJujo) 2 ' 

To understand what is going on let us consider a classical orbit with the center (X,Y). 
Then one can write 

x = X + v y /uj c , y = Y - v x /uj c . 

The quantity r c = v/uj c is the cyclotron radius of the orbit. We have two constants of 
motion, the energy E and X. In a long strip of width W the trajectories can be classified 
as a cyclotron orbits , skipping orbits, or traversing trajectory. In the (X, E) space such 
trajectories are separated by the line 



(X ± W/2f = r 



According to quantum mechanics, the motion is quantized, and we come to the following 
picture of quantum terms The cyclotron orbits (solid lines) correspond to Landau level, 
and they have zero group velocity. We have not shown traversing trajectories which cor- 
respond to higher energies. Dashed lines reproduce two sets of edge states (corresponding 
to skipping orbits). These sets carry the currents in opposite directions. 
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Figure 6.4: Typical electron trajectories in a 2D strip in magnetic field. 
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Figure 6.5: Electron terms in a 2D strip in magnetic field. 



If the Fermi level is situated between the Landau levels, only the edge states can con- 
tribute to the current. Their dispersion law can be obtained (approximately) from the 
Bohr-Sommerfeld quantization condition, 

h -1 (p p x dx + 7 = 27rn , n = 1, 2, ... 



One can show that for the rigid boundary the phase shift of the skipping orbit 7 = tt/2, 
while 

Px = mv x = (eH/c)(Y - y) . 

Thus, 

-4- £(Y -y)dx = 2n^- = 27i(n-—) . 

Consider an electron at the Fermi level. Its energy consists of {n — l/2)hu c (+ spin 
contribution which I do not discuss now), as well the the part 

E G = e F - (n - l/2)Hw c 
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due to electrostatic potential created by the edges, as well as by disorder. In an external 
potential, the center of the orbit moves with the drift velocity 

v d (R) = ^ [W(R) x H] 

which is parallel to the equipotentials. That can be easily shown from classical mechanics 
in a magnetic field. The typical spread of the wave function along the guiding center is 
spread within the range of magnetic length an- Now we can map the classical picture to 
quantum mechanics according to hk — > —x{eH/c). Thus, if the typical scale of the random 
potential of the disorder is greater than the magnetic length, we arrive at the picture. 



• Closed 

O Egde 






Wave vector, k «- 

Figure 6.6: Electron terms in the presence of long-range disorder. 



Assume that the edges are in a local equilibrium. Thus if there is a difference S( chemical 
potentials between the e<i<?es,then each channel contributes (e/h)5( to the current in the 
Hall direction. The system appears robust because to obtain a inter-channel exchange 
one needs tunneling with exponentially low probability. Actually we have an almost ideal 
ballistic conductor and the only difference with the systems discussed earlier is that the 
edge channels with different directions of the current do not overlap in space. 

In a typical realistic situation, the contacts are out of local equilibrium and the measured 
resistance depends on the properties of contacts. Consider for example, a situation when 
the edge channel at the lower edge are in equilibrium at chemical potential Ep, while the 
edge channel at the upper edge are not in local equilibrium. Then the current at the upper 
edge is not equipartitioned between N modes. Let f n is the fraction of the total current / 
that is carried by by states above Ep in the nth channel at the upper edge, I n = f n L The 
voltage contact at the upper edge will measure a chemical potential which depends on how 
it is coupled to each of the edge channels. The transmission probability T n is the fraction 
of the current I n that is transmitted through the voltage probe to a reservoir at chemical 
potential E F + 5(. The incoming current 



N 

hn = J2 T "f" I > With Z^» = 1 ' 

n n 



(6.2) 
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has to be balanced by an outgoing current, 



/out = UC{N - R) = UcY, T r 



(6.3) 



n 



since the voltage probe draws no net current. Thus the Hall resistance, 




) 



-i 



Rh — 



K h 



(6.4) 



el e 



n 



The Hall conductance remains quantized only if /„ = 1/N, or at T n = 1. The first case 
corresponds to local equilibrium, while the second case corresponds to an ideal contact. 
The Landauer-Buttiker formalism forms the basis on which anomalies in the QHE due 
to absence of local equilibrium in combination with non-ideal contacts can be treated 
theoretically. 

This is a simplified picture of the integer quantum Hall effect in the random potential. 
The real life is much more complicated. In particular, there exists an extremely interesting 
fractional quantum Hall effect, which manifests itself at fractional values of the filling 
factor. We do not discuss this effect in the present course. 

Role of localization 

As we have seen, at H = a 2D system with disorder should have its states localized at 
all energies. However, only extended states are sensitive to the flux and can provide the 
QHE. At the same time, ranges of energies with only localized states are needed to pin E F 
there and have finite plateaus. Thus, introduction of magnetic field must delocalize some 
states. As we have seen, extended modes appear near edges. However, extended states in 
a magnetic field are present also in the bulk of the conductor. 

To discuss this phenomenon let us recall the main relevant quantities. First, let us note 
that the condition 



for cyclotron motion is fully classical. In terms of quantum mechanical length, an = 



where N is the number of full Landau levels. The weak localization regime corresponds to 
the inequality 

a H < e, 

while the intermediate regime where a# <C r c while r c can be comparable with £ also exists. 



uj c t > 1, or r c < £, r c = v F /0 c 




the classical cyclotron radius r c can be written as 



kpa 2 H ~ an V 1 Epftw c ~ au^fN 



Strong magnetic field, uo c t ^> 1, r c I. 
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As we have discussed, in a uniform electric field the drift velocity directed along [E x H] 
appears, Vd = c(E/H). This concept can be generalized for the case of a smooth random 
potential V(x) with does not change a lot on the scale of cyclotron motion. Then, it r c is 
much less than the correlation length d of the random potential, the guiding center moves 
along the equipotential line V(r) = V. If its orbit is closed and embeds the area A(V), 
then the typical frequency is 



Ud __ 


U- 


1 _2c 


f dl dx± 


1 2c AV 


2tt ~" 


J v d _ 


= H 


[J dV _ 


H AA ' 



where dx± is an element of length in the direction of the potential gradient. Such a slow 
motion can be quantized for any Landau levels into locally equidistant levels with the 
separation hOd- The area between two quantized orbits is given by the relation 

HC 

HAA= — = $ ; AA = 2na 2 H . 

e 

Thus the flu of H in the area per state in a given Landau level corresponds to a flux 
quantum, as for free electron. 

Let us assume that the amplitude of the random potential is much less than hu c , so 
there is no inter-Landau-level mixing. Then the potential energy of a properly quantized 
levels should be added to ftw c (j + 1/2) for jth Landau band. The levels correspond to 
the orbits running around the potential "hills", or inside potential "lakes". Both kinds 
of states are localized. There is one and only one energy, E c , (in 2D case) at which 
the equipotential curves span the whole system (imagine filling up of the potential V(r) 
"terrain" with water). The characteristic size of the orbit, £ p , may be defined by the r.m.s. 
of the area enclosed by the equipotential contours for the localized states. It must blow up 
at E -> E e , 

& ~ \E - E c \- U * , v p >\. 

Such an orbit provides the way to transfer an electron between the edges. 
There is also a very interesting intermediate situation when 

aff -Cf <r c , or uj c t <ti 1 . 

As was shown by Khmelnitskii (1984), even in this region QHE plateaus can exist which 
are irrelevant to Landau levels. 

6.4 Fractional Quantum Hall Effect 

Fractional quantum Hall effect (FQHE) was discovered by the group of Tsui et Bell Lab- 
oratories [Hi]. Using a high-mobility GaAs/AlGaAs heterostructures they observed quan- 
tization of Hall conductance at filling factors v = 1/3 and 2/3 at very low temperatures 
(below 1 K). Later more rich structure, as shown in Figs. |6.7| and ^]8| at fractional fill- 
ing factors was discovered. It appears that only account of of Coulomb interaction leads 
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The fractional quan- 
tum Hall effect. The Hall resisti- 
vity (upper curve) of the in- 
version layer in a high mobili- 
ty AlGaAs/GaAs heterostructure 
shows plateaus at magnetic induc- 
tions B that correspond to the in- 
dicated filling factors v. At the 
same rilling factors, the magne- 
toresistivity (lower curve) shows 
minima 



100 150 200 250 BfkG 

Figure 6.7: 

to understanding of the problem. Now the studies of FQHE belong to the most active 
area of research. Below we shall provide a very brief sketch of the problems regarding 
electron-electron interaction in magnetic field and FQHE. 



Few electron with Coulomb interaction 

The role of electron-electron interaction is determined by the relation between the mean 
free distance between electrons, r s , and the Bohr radius, as = eh 2 /me 2 . At r s <C as one 
can use the usual mean field description of interacting electrons, considering screening, 
plasmons, charge density waves, etc. However, at r s > as the interaction energy becomes 
larger than the average kinetic energy. As a result, there exists a strong electron-electron 
correlation, and the electrons tend to crystallize. It is known that magnetic field enhances 
these effects. 

To get some understanding let us start with more simple problem of few electrons 
in a magnetic field. Historically, these studies appeared important because they led to 
discovery of a new state, the incompressible electron liquid, that is believed to transform 
into ( Wigner) crystal at very low densities. 
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MAGNETIC FIELD [T] 



Figure 6.8: Recent results on the fractional quantum Hall effect. 

Two electrons. Let us discuss the case of 2 electrons in a very strong magnetic field, 
^ c > e 2 /eau- This inequality means that Landau levels are not mixed by the Coulomb 
interaction (see below). 

Using symmetric gauge, A = (— Hy/2, Hx/2, 0)and introducing polar coordinates we 
easily obtain zeroth approximation Hamiltonian 

H = ^ V 2 | rnu 2 c p 2 | hw c d 



2m 8 2i dp> ' 

„ 2 d 2 Id Id 2 



dp 2 p dp p 2 dip 2 
This Hamiltonian commutes with the angular momentum, 

l _ k d 

i dip 

Thus it is natural to classify the states using the eigenvalues km of the angular momentum, 
l z . The eigenfunctions have the form 

*-M - |t> w ^ (~Q (4) ■ («) 
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Here n is non-negative integer, m is integer, Lj^' are Laguerre polynomials, while N nm = 
y / n!/27r2l m l(|m| + n) are normalization factors. The energy eigenvalues are 

E nm = huj c [n+(\m\-m + l)/2]. (6.6) 

The lowest Landau level corresponds to n — 0, m > 0. The Coulomb energy for the lowest 
state can be easily calculated as 

E c = (0m\e 2 /er\0m) = — + , 1/2) . (6.7) 

ea# m! 

At large m is decays as vrT 1 ! 2 . 

Two electrons are described by the Hamiltonian 

n (i) + n (2) + n int . 

It can be rewritten through the center-of-mass coordinate, R = (ri + r 2 )/v / 2, and the 
relative coordinate, r = (ri — r 2 )/v / 2, as 

n (R) + n (r) + n int (rV2). 

As a result, the center-of-mass motion is separated, and we look for a solution in the form 

*(R,r)=0(R)^(r). 

Now we are left with a single-particle problem for the relative motion. 

Since the interaction energy is radially-symmetric one can look for a solution in the 
form 

V>(r) = TZ(r)e~ imip 

with odd m because of the Pauli principle [ip(—r) = —ip(r)]. The radial Schrodinger 
equation is easily written in the dimensionless units as 

+ -[— + m+---)n = En, (6.8) 



2 dr 2 r dr 2 \r 2 4 r 

where r is measures in units of a#, E is measured in the units of hcu c , while dimensionless 
interaction constant is a = v / 2e 2 /ea^/ia; c . At large magnetic field this equation can be 
solved perturbatively with respect to a. In the lowest approximation we obtain: 

2 ea# m! 

The energy of the center-of-mass motion must be added. 

We find that the interaction destroys the degeneracy of the lowest Landau level. At 
large m the correction decreases because the electrons are less sensitive to interaction at 
long distances. 
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Three electrons. For 3 electrons we can also strip the center-of-mass motion. It can be 
done by the transform p = Or where 

/ 1/V2 -1/V2 
O = l/y/6 l/y/6 -2/y/6 
\ 1/V3 1A/3 1/V3 

After the transform the interaction Hamiltonian can be written as 

n int = ( - + 2 -y^- + ) . (6.10) 

eV2 \pi |Pi + V3p 2 | |v3p 2 -pi|/ 
Again, we can write the eigen function as a product 

*(Pi, P 2; Pa) = 0(p 3 MPi> P2) 

and in this way to reduce the problem to a two-particle one. 

An important point is that the probability density must be invariant under rotation 
about multiples of 7r/3. The resulting Hamiltonian also commutes with the total angular 
momentum, L. Then the states can be classified according to eigenvalues M of the orbital 
momentum. It was shown by R. Laughlin that a proper complete set to diagonalize a 
3-electron system can be written as 

\m,rri) = - [(z 2 + iz 1 ) 3m -(z 2 -iz 1 ) 3m ] (zf + z |)m' e -(M 2 +|2 2 | 2 ) _ ( 6 n ) 

Here Zi = ^ + T]i, eta are the Cartesian components of the vector p i /aH, F is a normal- 
ization factor. The states ( |6.11|) are the eigenstates of the total angular momentum with 



M = 3(m + m'). 

To diagonalize the system one has to solve the secular equation 

det \e5 mm 5 m > m > - (mm'\H int \mm')\ = . 

The crucial point is that the basis ( |6.11| ) is an extremely good starting approximation since 
off-diagonal elements of 7i int are typically at least 10 times less than the diagonal ones. 

The minimum angular momentum for which a non-degenerate solution exists is M = 
3 (to = 0, m' — 1). The next solution corresponds to M — 9, it is combined from the 
states (3,0) and (1,3). These states have the lowest energy at TC int = 0. The "charge 
density" for the state with M = 9 is 1/3 comparing to the state with M = 3. Since the 
angular momentum is conserved and the angular momentum corresponds to the area of an 
electronic state, the 3 electrons are "incompressible" . 



Fractional quantum Hall states 



It is impossible to diagonalize exactly the system of many electron states. An extremely 
effective approximate guess was suggested by R. Laughlin which we shall discuss for the 
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case of very large magnetic field when only the lowest Landau level is important. The 
single-electron states for that case can be written as 

( r | 0m ) = ^™ e -^ 2 / 4 , 
a H 

where z = x + iy. The complete set of iV-electron states with total angular momentum 
=1 m v are the Slater determinants 

N / JV \ 

..n)= £ (-i) p nv:^xphEw 2 ■ 

P(v 1 ...v N ) £t=l \ a=l / 

Since the ground state in the independent band approximation is a combination of Slater 
determinants, its general form is 

*(i . . . n) = n f( Zj - Zk ) ex P \ z «n ■ 

j<k \ a=l / 

There are several requirements to the functions f(z): 

• The function f(z) must be a polynomial of z; 

• Since $ should be a Fermion state, f(—z) = —f(z); 

• \I/ can be chosen as an eigenfunction of the total angular momentum. Therefore, the 
function f(z) has to be homogeneous. 

The simplest choice is 

f(z) = z m , (nodd). 
Thus the approximate wave function has the form 

*(i...iv)=n^-^) mex p(-iEi^i 2 ) • (° 2 ) 

j<k \ a=l / 

The Laughlin state ( |6.12|) describes a liquid-like system. The two-particle correlation func- 
tion 

r N 

(m) 
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J JJdr I/ |*(l...JV)| : 



u=3 

at small distances is proportional to \z± — Z2\ m that reflects the Pauli principle for the 
electrons. The smallest possible value of m is 3. The total angular momentum is just 
M = Nm, while the area covered by the electrons is A = N(2iTma 2 H ). Thus the average 
electron density is {2nma 2 H )~ 1 , and the filling factor is v = 1/m. To keep electrostatic 
stability one has to add the positive background. 
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The estimate for the Coulomb energy for the Laughlin state can me obtained as 

Because the correlation function decay strongly at small distances the incompressible liquid 
state appears more stable than the Wigner crystal. An interesting fact is that the Laughlin 
state appears the exact ground state for v — 1/m in the case of contact interaction, 
H int (r) (x5(r). 

As a consequence of the electron-hole symmetry it is easy to find the state corresponding 
to the filling factor (1 — v) if the state for v is known. 

Elementary excitations 

Elementary excitation are important both for transport and dynamics. Changing of energy 
of the electron system can be achieved by its compression, or, equivalently, by changing of 
angular momentum while keeping the neutralizing background. 

In other words, it means that new quasielectrons of quasiholes are introduced into the 
state ^ u if v ^ 1/m. An introduction of a quasihole can be represented as 

N 

^+ = A + (z Q )* u ( Zl ...z n ), A + {z ) = Y[{ Zj - z ) . 

i=i 

Let us estimate the effective charge of this excitation. The average area per particle which 
is covered in the state with filling factor v = 1/m is (N — \)(2nma 2 H ). It can be seen by 
direct calculation of the integral. The corresponding charge density is 

-Ne ^ e 

Po ~ (N- l){2-Kma 2 H ) ~ ~ 2ixma\ ' 

Thus, each electron occupies the area with m flux quanta. Its charge must be compensated 
by a positive background. 

In the state the maximum angular momentum per particle is increased by 2-Ka 2 H . 
This corresponds to the change in the charge density which is equivalent to the positive 
charge +e/m. 

Quasielectrons can be created in a similar way, 

N ( d 

% = A-(z Q )m u { Zl ...z n ), A + (z ) = Yl - z { 

Here the partial derivative acts only on the polynomial part of the wave function ^ u leaving 
alone the Gaussian part. It can be shown that the effective charge of the quasielectron is 
—e/m. 

The gaps between the ground and excited states were observed directly from tempera- 
ture dependences of conductance. It appears that the quasiparticles can be considered as 
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particle with so-called fractional statistics - anyons. Very interesting collective excitations 
were also predicted and observed experimentally by inelastic light scattering. H. Stormer, 
D. Tsui and R. Laughlin were awarded by the Nobel Prize 1998 for their discovery of 
FQHE. 

However, the story is not over. Very specific features of Hall conductance were observed 
at v — p/q where p, q are integers, both for odd and even denominator q. These features 
were not explained by original theory by Laughlin. It appears, that at odd denominators 
the electrons also condense in some quantum liquids. However, the properties of that 
liquids differ significantly from those of the incompressible Laughlin liquid. 

The above discussion is definitely not an "explanation" of FQHE. It just demonstrates 
some basic trends in the field. More work has to be done to understand the whole physical 
picture and to construct a proper transport theory. 



Chapter 7 

Noise in mesoscopic systems 



Shot noise, the time-dependent fluctuations in the electrical current due to the discreteness 
of the electron charge. 

The shot-noise power in small conductors possesses special features which are important 
both for understanding correlations phenomena in mesoscopic systems and for application 
of mesoscopic devices. The noise is maximal if the electron transmission through the system 
is fully uncorrelated. This maximal value is called the Poisson limit. It mesoscopic systems 
it can be suppressed as a result of correlations in the electron transmission imposed by the 
Pauli principle. This suppression takes on simple universal values in a symmetric double- 
barrier junction (suppression factor |), a disordered metal (factor |), and a chaotic cavity 
(factor |). Loss of phase coherence has no effect on this shot-noise suppression, while 
thermalization of the electrons due to electron-electron scattering increases the shot noise 
slightly. Sub-Poissonian shot noise has been observed experimentally. So far unobserved 
phenomena involve the interplay of shot noise with the Aharonov-Bohm effect, Andreev 
reflection, and the fractional quantum Hall effect. 

Below we discuss these processes following the review Ref. p3 . 



7.1 Current fluctuations 

In 1918 Schottky reported that in ideal vacuum tubes, where all sources of spurious noise 
had been eliminated, there remained two types of noise in the electrical current, described 
by him as the Warmeeffekt and the Schroteffekt. The first type of noise became known 
as Johnson- Nyquist noise, or simply thermal noise. We discussed this type of noise in 
connection with fluctuation-dissipation theorem. It is due to the thermal motion of the 
electrons and occurs in any conductor. The second type of noise is called shot noise, caused 
by the discreteness of the charge of the carriers of the electrical current. Not all conductors 
exhibit shot noise. 

Noise is characterized by its spectral density or power spectrum P(oo), which is the 
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Fourier transform at frequency uj of the current-current correlation function, 

oo 

P(uj) = 2 J dte iuJt (AI(t + t )AI(t )) . (7.1) 

— oo 

Here AJ(t) denotes the time- dependent fluctuations in the current at a given voltage V 
and temperature T. The brackets (• • • ) indicate an ensemble average or, equivalently, an 
average over the initial time to- Both thermal and shot noise have a white power spectrum 
- that is, the noise power does not depend on uj over a very wide frequency range. Shot 
noise (V ^ 0, T = 0) is more interesting then the thermal one, because it gives information 
on the temporal correlation of the electrons, which is not contained in the conductance. 
In devices such as tunnel junctions, Schottky barrier diodes, p-n junctions, and thermionic 
vacuum diodes, the electrons are transmitted randomly and independently of each other. 
The transfer of electrons can be described by Poisson statistics, which is used to analyze 
events that are uncorrelated in time. For these devices the shot noise has its maximum 
value 

P = 2el = P Poisson , (7.2) 

proportional to the time-averaged current /. Q Correlations suppress the low-frequency 
shot noise below Pp isson- One source of correlations, operative even for non-interacting 
electrons, is the Pauli principle, which forbids multiple occupancy of the same single- 
particle state. A typical example is an ideal ballistic point contact, where P = because 
the stream of electrons is completely correlated by the Pauli principle in the absence of 
impurity scattering. 

Progress in nanofabrication technology has revived the interest in shot noise, because 
nanostructures allow measurements to be made on "mesoscopic" length scales that were 
previously inaccessible. The mesoscopic length scale is much greater than atomic dimen- 
sions, but small compared to the scattering lengths associated with various inelastic pro- 
cesses. Mesoscopic systems have been studied extensively through their conductance. Noise 
measurements are much more difficult, but the sensitivity of the experiments has made a 
remarkable progress in the last years. 

To analyze the noise properties we shall apply scattering formalism. Incoming and 
outgoing waves are specified as shown in Fig. [771] Each lead contains N incoming and N 
outgoing modes at energy The incoming and outgoing modes are related by a 2N x 2iV 
scattering matrix S 




where Ii,Oi,l2, O2 are the iV-component vectors denoting the amplitudes of the incoming 
(J) and outgoing (O) modes in lead 1 and lead 2. The scattering matrix can be decomposed 

Equation ( |7.2| ) is valid for uj < t -1 , with r the width of a one-electron current pulse. For higher 
frequencies the shot noise vanishes. 

2 We assume only elastic scattering so that energy is conserved. 
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Figure 7.1: Schematic representation of the transport through the conductor. Incoming 
states (/) are scattered into outgoing states (O), by a scattering region (dashed). A cross 
section in lead 1 and its coordinates are indicated. 



in N x N reflection and transmission matrices 



Sn s 12 \ _ / r t' 
S21 s 22 j " I t r' 



(7.4) 



where the N x N matrix Sb a contains the amplitudes Sb n ,am from incoming mode m in 
lead a to outgoing mode n in lead b. Because of flux conservation S is a unitary matrix. 
Moreover, in the presence of time-reversal symmetry S is symmetric. 
The current operator in lead 1 is given by 



oo oe> 



/(*) = J d£ J d£> J ^ £ > £ ') ^(e) a p {e') e u ^ H , (7.5) 

<*,0 

where a) a (e) [a a {e)\ is the creation [annihilation] operator of scattering state ip a (r,e). We 
have introduced the indices a = (a,m), (3 = (b, n) and the coordinate r = (x,y). The 
matrix element I a p(e,e') is determined by the value of the current at cross section Si in 
lead 1, 

Iap(e,e') = \ Jdy{ip a (r,e)[v x ip p (r,e')]* + ip* p (r,e')v x ip a (r,e)} . (7.6) 

Si 

Here, v x is the velocity operator in the x-direction. At equal energies, Eq. ( |7.6|) simplifies 
to 

N 

(e) . (7.7) 

The average current follows from 

(al(e)d (e'))=8 a p8(e-e')f a (e), (7.8) 

where f a is the Fermi-Dirac distribution function in reservoir a: 

fi(e) = f(e — Ep — eV) , (9a) 
f 2 (e) = f(e-E F ), (9b) 
f{x) = [l + exp^/A^T)]- 1 , (9c) 
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with Fermi energy Ep. The result is f] 

('(*)> = r > ] ldef a (e)I aa (e,e) = Ids [f x {e) - / 2 (e)] Trt^t^) , (7.10) 



h 



Jdef a (e)I aa (e,e) = & - Jde [f x {e) - f 2 (e)] T*t(e)tt(e) 



where we have substituted Eq. (|7.7|) and used the unitarity of S. The linear-response 
conductance, G = \im v ^o(I) /V , becomes 



oo 

,2 



G= 6 - [ de(- d 4-)Trt{e)t\e). (7.11) 



h J \ ds 
o 

At zero temperature we reproduce the Landauer formula (per spin) 

2 2 N 

G=^Trttt = ^T n . (7.12) 

n=l 

Here t is taken at Ep and T n e [0, 1] is an eigenvalue of t^. The conductance is thus fully 
determined by the transmission eigenvalues. Knowledge of the transmission eigenstates, 
each of which can be a complicated superposition of incoming modes, is not required. 

In order to evaluate the shot-noise power we substitute the current operator (|7.5|) into 
Eq. ( |7.1|) and determine the expectation value. It can be shown from direct quantum- 
mechanical calculation that 

(a\a 2 ala 4 ) - {a\a 2 ) {a\a^} = 5i4<Wi(l - $2) = A1234 , (7.13) 

where e.g. 5 12 stands for 8 a p8(e— e'). Equation ( |7. 13| ) shows that there are cross correlations 
between different scattering states. Although this bears no effect on the time-averaged 
current, it is essential for the current fluctuations. For the noise power at low-frequancy 
limit one finds, 

00 

P = 2^|&{[/ 1 (l-/ 2 ) + / 2 (l-/ 1 )]Trttt(l-ttt) 



+ [A(l - /1) + /a(l - /a)] Tttttttt} , (7.14) 

where we have again used the unitarity of S. 

Equation (|7.14|) allows us to evaluate the noise for various cases. Below we will assume 
that eV and fc^T are small enough to neglect the energy dependence of the transmission 
matrix, so that we can take t at e — Ep. Let us first determine the noise in equilibrium, 
i.e. for V = 0. Using the relation f(l — f) = —ksTdf/de we find 

e 2 e 2 A 

P = 4A; B T-Trtt t = 4A; B T- ^T n , (7.15) 

71=1 



3 Note that in the above derivations the absence of spin and valley degeneracy has been assumed for 
notational convenience. It can be easily included. 



7.2. TWO SIMPLE APPLICATIONS 
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which is indeed the Johnson- Nyquist formula. For the shot-noise power at zero temperature 
we obtain 



,2 



A' 



P = 2eV^Trtt t (l-tt t ) =2eV^ £ T„(l - T n ) . (7.16) 

n=l 

One notes, that P is again only a function of the transmission eigenvalues. 

It is clear from Eq. ([7.161) that a transmission eigenstate for which T n = 1 does not 
contribute to the shot noise. This is easily understood: At zero temperature there is a non- 
fluctuating incoming electron stream. If there is complete transmission, the transmitted 
electron stream will be noise free, too. If T n decreases, the transmitted electron stream 
deviates in time from the average current. The resulting shot noise P is still smaller than 
-Ppoisson, because the transmitted electrons are correlated due to the Pauli principle. Only 
if T n <C 1, the transmitted electrons are uncorrelated, yielding full Poisson noise. 

The generalization of Eq. ( [7.1 6|) to the non-zero voltage, non-zero temperature case is 

2 N 

P = 2^-Y1 [ 2 ^ TT « + T «( X - T n )eV coth{eV/2k B T)} . (7.17) 



n=l 



The crossover from the thermal noise (|7.15| ) to the shot noise ( [7.1 6|) depends on the trans- 
mission eigenvalues. 



7.2 Two Simple Applications 

The above results are valid for conductors with arbitrary (elastic) scattering. If the trans- 
mission eigenvalues are known, the conduction and noise properties can be readily calcu- 
lated. Below, this is illustrated for two simple systems. More complicated conductors are 
discussed in Sees. OHZ2H. 



Tunnel barrier 

In a tunnel barrier, electrons have a very small probability of being transmitted. We model 
this by taking T n -C 1, for all n. Substitution into the formula for the shot noise (|7.16|) and 
the Landauer formula for the conductance ( [7.1 2| ) yields P = Pp isson at zero temperature. 
For arbitrary temperature we obtain from Eq. ( |7.17| ), 

P = coth(eV/2A; B T) P Poisson . (7.18) 

This equation, describes the crossover from thermal noise to full Poisson noise. For tunnel 
barriers this crossover is governed entirely by the ratio eV/ksT and not by details of the 
conductor. This behavior has been observed in various systems. 
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Figure 7.2: (a) Conductance G (dashed line) and shot-noise power P (full line) versus 
Fermi energy of a two-dimensional quantum point contact, according to the saddle-point 
model, with uj y = 3u x . (b) Experimentally observed G and P versus gate voltage V g 
(unpublished data from Reznikov et al. similar to the experiment of Ref. [B3], but at a 
lower temperature T = 0.4 K). 



Quantum Point Contact 

As we know, the conductance displays a stepwise increase in units of Go as a function of V g . 
At the conductance plateaus the shot noise is absent, as follows from Eq. (|7.16|) . However, 
in between the plateaus, where the conductance increases by Go, there is a transmission 
eigenvalue which is between and 1. As a consequence, the shot noise has a peak. Results 
for the conductance and the shot-noise power are displayed in Fig. |7T^| a. The shot noise 
peaks in between the conductance plateaus and is absent on the plateaus. For large N, the 
peaks in the shot noise become negligible with respect to the Poisson noise, in agreement 
with the classical result. 

The prediction of this quantum size-effect in the shot noise formed a challenge for 
experimentalists. Recent experiments at high frequencies by Reznikov et al. [<VZ] have 
unambiguously demonstrated the occurrence of suppressed shot noise on the conductance 
plateaus. Experimental data of Reznikov et al. are shown in Fig. [7.2| b. 



7.3 Phase Breaking, Thermalization, 
and Inelastic Scattering 

Noise measurements require rather high currents, which enhance the rate of scattering 
processes other than purely elastic scattering. The phase-coherent transmission approach 
is then no longer valid. Below, we discuss a model in which the conductor is divided 
in separate, phase-coherent parts connected by charge-conserving reservoirs. This model 
includes the following types of scattering: 

• Quasi- elastic scattering. Due to weak coupling with external degrees of freedom the 
electron-wave function gets dephased, but its energy is conserved. In metals, this 
scattering is caused by fluctuations in the electromagnetic field. 
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• Electron heating. Electron-electron scattering exchanges energy between the elec- 
trons, but the total energy of the electron gas is conserved. The distribution func- 
tion is therefore assumed to be a Fermi-Dirac distribution at a temperature above 
the lattice temperature. 

• Inelastic scattering. Due to electron-phonon interactions the electrons exchange en- 
ergy with the lattice. The electrons emerging from the reservoir are distributed 
according to the Fermi-Dirac distribution, at the lattice temperature T. 




f 12 (e) 



Figure 7.3: Additional scattering inside the conductor is modeled by dividing it in two 
parts and connecting them through another reservoir. The electron distributions in the left 
and the right reservoir, /i(e) and f 2 (s), are Fermi-Dirac distributions. The distribution 
/12(e) in the intermediate reservoir depends on the type of scattering. 

The model is depicted in Fig. [7.3| . The conductors 1 and 2 are connected via a reservoir with 
distribution function /12(e)- The time-averaged current I m through conductor m = 1,2 is 
given by 

h = (Gt/e) Jde[h{e)-h2{e)], (19a) 
J 2 = (G 2 /e) jde[f 12 (e)-f 2 (e)}. (19b) 

The conductance G m = 1/R m = Gq ^2^=i Tn m \ with T^ 1 the n-th transmission eigenvalue 
of conductor m. We assume small eV and ksT, so that the energy dependence of the 
transmission eigenvalues can be neglected. 

Current conservation requires that I\ — I 2 = I. The total resistance of the conductor 
is given by Ohm's law, 

R = R 1 + R 2l (7.20) 

for all three types of scattering that we consider. 

The time- averaged current (|7.19|) depends on the average distribution /12(e) in the 
reservoir between conductors 1 and 2. In order to calculate the current fluctuations, we 
need to take into account that this distribution varies in time. We denote the time- 
dependent distribution by fi 2 (e,t). The fluctuating current through conductor 1 or 2 
causes electrostatic potential fluctuations 8(j)i 2 (t) in the reservoir, which enforce charge 
neutrality. Assuming that the reservoir has a Fermi-Dirac distribution f\ 2 (e,t) = f[e — 

4 The model is not suitable for transport in the ballistic regime or in the quantum Hall effect regime, 
where a different type of "one-way" reservoirs is required. 
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Ep — eVw — e5<pi2{t)], with Ep + eVu the average electrochemical potential in the reservoir. 
As a result, it is found that the shot-noise power P of the entire conductor is given by 

R 2 P = i^Pj + R\P 2 . (7.21) 

In other words, the voltage fluctuations add. The noise powers of the two segments depend 
solely on the time-averaged distributions, 

P m = 2G m Jde [f m (l - f m ) + Ml - /12)] + 2S m Jde (f m - f l2 ) 2 , (7.22) 

where S m = Go 2^n=i Tn m \l — T„ ). The analysis first given in Ref. [ESI] is easily gener- 
alized to arbitrary distribution / 12 . Then, we have fi 2 (e,t) = fi 2 [e — e5(fii 2 (t)}. It follows 
that Eqs. ( |7.21|) and (|7.22|) remain valid, but j 12 (^) may be different. Let us determine the 
shot noise for the three types of scattering. 

Quasi-elastic scattering. Here, it is not just the total current which must be conserved, 
but the current in each energy range. This requires 

, ( ,_ G 1 f 1 (e) + G 2 f 2 (e) 

Me) - . (7.23) 

We note that Eq. ( ^23|) implies the validity of Eq. ( |7720|) . Substitution of Eq. (|?T23|) into 
Eqs. ( [7.2 1|) and ( |7.22|) yields at zero temperature the result: 

P = (R1S1 + R 2 S 2 + RiR\ + R\R 2 ) R~ 3 -Ppoisson • (7.24) 

Electron heating. We model electron-electron scattering, where energy can be exchanged 
between the electrons at constant total energy. We assume that the exchange of energies 
establishes a Fermi-Dirac distribution fi 2 {s) at an electrochemical potential Ep + eV\ 2 and 
an elevated temperature T 12 . From current conservation it follows that 

V l2 = (R 2 /R) V . (7.25) 

Conservation of the energy of the electron gas requires that T 12 is such that no energy is 
absorbed or emitted by the reservoir. This implies 

m o m o Rt Ro ,_ . 

T? 2 = T 2 + — , (7.26) 

with the Lorentz number Co = ^(nks/ 'e) 2 '. At zero temperature in the left and right 
reservoir and for R l = R 2 we have k B T 12 = (V3/2n)eV ~ 0.28eV. For the shot noise at 
T = 0, we thus obtain using Eqs. (|7.21|) and (|7.22|) the result: 

R1S1 + R 2 S 2 + ^y/3RiR 2 R\(l — R1S1) + R 2 (l — R 2 S 2 ) 
+ 2R\S 1 In (l + e-^V^i/ 3 ^ 

+ 2R\S 2 In (l + e-V^T^)] | R -2 p poisson (7 27) 
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Inelastic scattering. The distribution function of the intermediate reservoir is the Fermi- 
Dirac distribution at the lattice temperature T, with an electrochemical potential /ii 2 = 



Ep + eVi2i where Via is given by Eq. ( |7.25|) . This reservoir absorbs energy, in contrast to 



the previous two cases. The zero-temperature shot-noise power is given by 

P = (R 3 ^ + R 3 2 S 2 ) R- 2 P Poisson . (7.28) 

This model will be applied to double-barrier junctions, and disordered conductors in the 
following sections. Quite generally, we will find that quasi-elastic scattering has no effect 
on the shot noise, while electron heating leads to a small enhancement of the shot noise. 
Inelastic scattering suppresses the shot noise in most cases, but not in the double-barrier 
junction. 



7.4 Double-Barrier Junction 

Resonant Tunneling 

Below we will only consider the zero-frequency, low- voltage limit, in order to treat the 
double-barrier junction on the same footing as the other systems described in these lectures. 
We assume high tunnel barriers with mode-independent transmission probabilities ri, r 2 <C 
1. 

The transmission eigenvalues through the two barriers in series, as we have derived 
above, can be re-written as, 

T n = 2 , (7.29) 

2-r 1 -r 2 -2 v / i-r 1 -r 2 cos0 n v ' 

where <p n is the phase accumulated in one round trip between the barriers. The density 
p(T) = (Yin ${T—T n )) of the transmission eigenvalues follows from the uniform distribution 
of 4> n between and 2tt, 

p(T) = 1 ' , = , Te[T_,T + ], (7.30) 

HK ' 7r(ri + r 2 ) V /T 3 (T+ - T) y ' 



p{T) = otherwise, with T_ = rir 2 /vr 2 and T+ = 4r 1 r 2 /(ri + T 2 ) 2 . The density ( fT\30[) is 
plotted in Fig. [7T% . 

The average conductance, 

l 

(G) = G J dTp(T) T = G N^^r , (7.31) 
o 

is just the series conductance of the two tunnel conductances. The resonances are averaged 
out by taking a uniform distribution of the phase shifts <p n . Physically, this averaging 
corresponds either to an average over weak disorder in the region between the barriers, 
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Figure 7.4: The distribution p(T) of transmission eigenvalues T for (a) a double-barrier 
junction, according to Eq. ( |7.29|) with I\ = T 2 = 0.01; (b) a chaotic cavity, which we do 
not discuss here in detail; and (c) a disordered wire, according to Eq. (|7.35|) with L = 20£. 
Each structure has a bimodal distribution. 



or to a summation over a large number of modes if the separation between the barriers is 
large compared to the Fermi wave length, or to an applied voltage larger than the width 
of the resonance. For the shot-noise power one obtains 

(P) = P J dT p(T) T(l — T) = ^l±lLp PoissoQ , (7.32) 
o 

using Eqs. (|7.30|) and (|7.31|) . For asymmetric junctions, one barrier dominates the transport 
and the shot noise equals the Poisson noise. For symmetric junctions, the shot noise gets 



suppressed down to (P) = ^Pp isson for Ti = T 2 . The theoretical result (7.32) is in 
agreement with the several experimental observations. 

The suppression of the shot noise below -Pp i SS on in symmetric junctions is a consequence 
of the bimodal distribution of transmission eigenvalues, as plotted in Fig. |7.4| a. Instead of 
all T n 's being close to the average transmission probability, the T ra 's are either close to 
or to 1. This reduces the sum T n (l — T n ). A similar suppression mechanism exists for shot 
noise in chaotic cavities and in disordered conductors. 

Phase coherence is not essential for the occurrence of suppressed shot noise. The method 
described above (with G m = S m = G NT m for m = 1,2) shows that both quasi-elastic 



scattering [see Eq. ( |7.24|) 1 and inelastic scattering [see Eq. ( |7.28|) 1 do not modify Eq. (|7.32|) . 
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Thermalization of the electrons in the region between the barriers enhances the shot noise, 
as follows from Eq. (|7.27 ). For I\ = T2 we find 



P 



+ ^ In (l + 



-^Poisson — 0.58-Ppoisson > 



(7.33) 



which is slightly above the one-half suppression in the absence of thermalization. 



Coulomb Blockade 

The suppression of the shot noise described in the previous Section is due to correlations 
induced by the Pauli principle. Coulomb interactions are another source of correlations 
among the electrons. As we have discussed, a measure of the importance of Coulomb 
repulsion is the charging energy Ec = e 2 /2C of a single electron inside the conductor with 
a capacitance C. If eV < E c , conduction through the junction is suppressed. At eV > E c , 
one electron at a time can tunnel into the junction. The next electron can follow, only after 
the first electron has tunneled out of the junction. This is the single-electron tunneling 
regime. 

Experiments on the noise suppression in Coulomb-blockade regime were reported in 
Reef. [03]. Here, the double-barrier junction was formed by a scanning-tunneling micro- 
scope positioned above a metal nanoparticle on an oxidized substrate. Due to the small 
size of the particle, E c > lOOO/c^T, at T = 4 K. The relative heights of the two tunnel 
barriers can be modified by changing the tip-particle distance. Experimental results for 
an asymmetric junction are plotted in Fig. |7.5| . The I-V characteristics display a stepwise 
increase of the current with the voltage. (Rotating the plot 90° yields the usual presenta- 
tion of the 'Coulomb staircase.') At small voltage, / ~ due to to the Coulomb blockade. 
At each subsequent step in /, the number of excess electrons in the junction increases by 
one. The measured shot noise oscillates along with the step structure in the I-V curve. 
The full shot-noise level P = Pp isson is reached at each plateau of constant /. In between, 
P is suppressed down to |Pp i S son- The experimental data are in excellent agreement with 
the theory. 

A qualitative understanding of the periodic shot-noise suppression caused by the Coulomb 
blockade goes as follows: On a current plateau in the I-V curve, the number of electrons in 
the junction is constant for most of the time. Only during a very short instance an excess 
electron occupies the junction, leading to the transfer of one electron. This fast transfer 
process is dominated by the highest tunnel barrier. Since the junction is asymmetric, Pois- 
son noise is expected. The situation is different for voltages where there is a step in the I-V 
curve. Here, two charge states are degenerate in total energy. If an electron tunnels into 
the junction, it may stay for a longer time, during which tunneling of the next electron is 
forbidden. Both barriers are thus alternately blocked. This leads to a correlated current, 
yielding a suppression of the shot noise. 

An essential requirement for the Coulomb blockade is that G < e 2 /h. For larger G the 
quantum-mechanical charge fluctuations in the junction become big enough to overcome 
the Coulomb blockade. The next Section will discuss shot noise in a quantum dot, without 
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Figure 7.5: Experimental results by Birk et al. [ED] in the single-electron tunneling regime. 
The double-barrier junction consists of a tip positioned above a nanoparticle on a substrate, 
(a) Experimental voltage V versus current /. (b) Shot-noise power P versus I. Squares: 
experiment; solid line: theory. 

including Coulomb interactions. This is justified as long as G > e 2 /h. For smaller G, the 
quantum dot behaves essentially as the double-barrier junction considered above. 

Disordered Metal 
One-third suppression 

We now turn to transport through a diffusive conductor of length L much greater than 
the mean free path £, in the metallic regime (L <C localization length). The average 
conductance is given by the Drude formula, 

N£ 

(G) = G — , (7.34) 

up to small corrections of order Go (due to weak localization). The mean free path £ = £tr 
equals the transport mean free path £ tr times a numerical coefficient, which depends on 
the dimensionality d of the Fermi surface (02 = vr/2, a 3 = 4/3). 

From Eq. ( |7.34j) one might surmise that for a diffusive conductor all the transmission 
eigenvalues are of order £/L, and hence 1. This would imply the shot-noise power P = 
-fpoisson of a Poisson process. This surmise is completely incorrect. A fraction £/L of the 
transmission eigenvalues is of order unity (open channels), the others being exponentially 
small (closed channels). For f<L<C N£, the density of the T n 's is given by 




(7.35) 
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p(T) = otherwise, with T_ = 4e 2L / e . The density p(T), plotted in Fig. 7_Ac, is again 
bimodal with peaks near unit and zero transmission. 



One easily checks that the bimodal distribution (|7.35 ) leads to the Drude conductance 
([7.34] ) . For the average shot-noise power it implies 



(P) = Po 



N£ 
3L 



Poisson 



(7.36) 



This suppression of the shot noise by a factor one-third is universal, in the sense that it 
does not depend on the specific geometry nor on any intrinsic material parameter (such as 



Dependence on wire length 



The one-third suppression of the shot noise breaks down if the conductor becomes too short 
or too long. Upon decreasing the length of the conductor, when L becomes comparable to 
£, the electron transport is no longer diffusive, but enters the ballistic regime. Then the 
shot noise is suppressed more strongly, 



P = Ml-(l + L/£)- 3 ]P ] 



Poisson 



(7.37) 



For L <C I there is no shot noise, as in a ballistic point contact. Equation ( [7.37 ) is exact 
for a special model of one-dimensional scattering, but holds more generally within a few 
percent. The crossover of the shot noise from the ballistic to the diffusive regime is plotted 
in Fig. |7.6| . Upon increasing L at constant cross section of the conductor, one enters the 
localized regime. Here, even the largest transmission eigenvalue is exponentially small, so 
that P = Ppoisson- Experimentally, the crossover from the metallic to the localized regime 
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Figure 7.6: The shot-noise power P of a disordered metallic wire as a function of its 
length L, as predicted by theory. Indicated are the elastic mean free path £, the electron- 



electron scattering length l e 
are interpolations. 



and the electron-phonon scattering length l ep . Dotted lines 



is usually not reached, because phase coherence is broken when L is still much smaller 
than the localization length N£. In the remainder of this Section, we apply the method of 
Sec. [7TB] to determine the effect of phase breaking and other inelastic scattering events on 
the shot noise in a disordered metal. We divide the conductor into M segments connected 
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by reservoirs, taking the continuum limit M — > oo. The electron distribution at position 
x is denoted by f(e, x). At the ends of the conductor f(e, 0) = fi(e) and f(e, L) = /2(e), 
i.e. the electrons are Fermi-Dirac distributed at temperature T and with electrochemical 
potential /x(0) = E F + eV and fi(L) = E F , respectively. It follows from Eqs. (|7.21|) and 
( |7.22| ) that the noise power is given by 



4 
RL 



dx / de f(e, x) [1 — f(e, x)] . 



(7.38) 



We evaluate Eq. ( |7.38|) for the three types of scattering discussed in Sec. [7]3 



Quasi-elastic scattering. Current conservation and the absence of inelastic scattering 
requires 

f(e,x) = I ^f(e,0) + yf(e,L). (7.39) 



The electron distribution at x = L/2 is plotted in the inset of Fig. [7/7|. Substitution of Eq. 
(fT39|) into Eq. Q7gg) yields 



[Ak B TG + e/coth(eV/2A; B T)] 



(7.40) 



At zero temperature the shot noise is one-third of the Poisson noise. The same result 
follows from the phase-coherent theory [Eqs. (|7.17|) and ( 7.35 )], demonstrating that quasi- 
elastic scattering has no effect on the shot noise. The temperature dependence of P is 
plotted in Fig. [7T7| . 

Electron heating. The electron-distribution function is a Fermi-Dirac distribution with 
a spatially dependent electrochemical potential /i(x) and temperature T e (x), 



fj,(x) 
TJx) 



1 + exp 
L - 



e — n(x) 
k B T e (x) 



E F + 



x 



-eV 



(41a) 

(41b) 
(41c) 



= y/T* + (x/L)[l-(x/L)]VyCo, 
cf. Eqs. (|7.25|) and (|7.26|) . Equations (|7.38|) and ( |7.41| ) yield for the noise power the result 



P = 2k B TG + 2el 



2tt fk B T 

7! 



+ 



2tt 



arctan 



V3 eV 



2tt k B T 



plotted in Fig. |7.7| . In the limit eV ^> k B T one finds 

P = !v3Ppoisson — 0.43 Ppoisson • 



(7.42) 



(7.43) 



Electron-electron scattering increases the shot noise above |Pp i SS on because the exchange 
of energies makes the current less correlated. 
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Figure 7.7: The noise power P versus voltage V for a disordered wire in the presence 
of quasi-elastic scattering [solid curve, from Eq. ( |7.40|) 1 and of electron heating [dashed 
curve, from Eq. ( |7.42|) 1. The inset gives the electron distribution in the middle of the 
wire at A;#T = ^eV. The distribution for inelastic scattering is included for comparison 
(dash-dotted). Experimental data of Steinbach, Martinis, and Devoret [EEI] on silver wires 
at T = 50 mK are indicated for length L = lfim (circles) and L = 30/im (dots). 



Inelastic scattering. The electron-distribution function is given by 



f(e,x) 



1 + exp 



e — fi(x) 



(7.44) 



with n(x) according to Eq. (|41b| ). We obtain from Eqs. ( |7.38| ) and (|7.44|) that the noise 
power is equal to the Johnson- Nyquist noise for arbitrary V. The shot noise is thus 
completely suppressed by inelastic scattering. 

The dependence of the shot-noise power on the length of a disordered conductor is 
plotted in Fig. |7.6| . The phase coherence length (between £ and l ee ) does not play a role. 

Results of accurate experiments by Steinbach, Martinis, and Devoret [ESI] on silver wires 
are shown in Fig. \T7\. The noise in a wire of L = 30 /im is in excellent agreement with 
the hot-electron result (|7.42|) . For the L = 1 /im wire the noise crosses over to the elastic 
result ( |7.40| ), without quite reaching it. 
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Appendix A 
Second quantization 



Now we briefly discuss the way to describe many-electron states with the occupation num- 
bers. 



Bosons 

In general case, the total wave function of bosons is symmetric in replacement of the 
particles. Thus it can be expressed as a symmetric product of individual wave functions 

{ N IN ! \ 

$jviiv 2 ... = y 1 tvT y £^w(£i)v»(&)---Vp*(6v), ( A - x ) 

where Pi label the states, ip Pi are while the sum is calculated over all the permutations of 
{pi}. The coefficient is introduced to make the to total function normalized: J |$| 2 Yli d£,i = 
1. 

The main idea is to consider $ jv function of the occupation numbers Aj. Assume 

that we have an arbitrary one-particle symmetric operator 

^ (1) = E/i 1} ( A - 2 ) 

a 

where f a acts only upon the functions of £ a . In is clear that acting upon ^n 1 n 2 ... h changes 
the state of only one particle. So in is reasonable to introduce the operators with matrix 
elements 

(bi)^ 1 = VNi , (bl)^ = [ib^ 1 ] * = (A.3) 

It is just annihilation and creation operators introduced earlier. The operator ( A.2|) can 
be expressed in terms of the creation-annihilation operators as 

F {1) = Y,f^ ( A - 4 ) 

ik 

where 

/i 1} = / ^(0/ (1 V(0^- (A.5) 
109 
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One can easily prove this relation comparing both diagonal and off-diagonal matrix ele- 
ments of the operators. A 2-particle symmetric operator 



^ (2) = £/- (A.6) 



a.b 



(2) 

where fv acts upon the functions of the variables £ a and can be expressed as 



^ (2) = £/K^An (A.7) 



iklm 



where 



Electrons 



ft = / ^(6K(&)/ (2 VKa)M6) dCi d£ a . (A.8) 



Now we turn to the Fermi statistics to describe electrons. According to the Pauli principle, 
the total wave function should be anti- symmetric over all the variables. So the occupation 
numbers could be or 1. In this case we get instead of (|A.1|) 



^N lN2 ... = -7= £(-l) P ^ Pl (ei)^ P2 (6) • • • <PpA&) ( A - 9 ) 

ViV. p 

where all the numbers pi, P2, ■ ■ ■ , Pn are different. The symbol (— l) p shows that odd and 
even permutations enter the expression ( |A.9|) with opposite signs (we take the sign '+' for 
the term with pi < p 2 < ... < Pn)- Note that the expression ( |A.9|) can be expressed as the 
determinant of the matrix with the elements = (1/ 'y/N\)<p Pi (£k) which is often called 
the Slater determinant. 

The diagonal matrix elements of the operator are 

^ (1) = E^ ( A - 1Q ) 

i 

just as for the Bose particles. But off-diagonal elements are 

where the sign is determined by the parity of the total number of particles in the states 
between the i and k ones. Q Consequently, for Fermi particles it is convenient to introduce 
the annihilation and creation operators as 



!)J = (-1)E&* (A.12) 



1 Note that for Bose particles similar matrix elements are (F^J ^.'f^^ — fik^ V^i^k 



Ill 

We immediately get (Check!) the commutation rules for the Fermi operators: 

1 = a { a\ + a\a k = 8 i} 



Hk 

„t„t ' 



{ ai a k } = |aX|=0. (A.13) 

The product of Fermi operators are 

a\ai = Ni , aioj = 1 - Ni . (A. 14) 

One can express all the operators for Fermi particles exactly in the same way as the Bose 
ones, Eqs. flOR , (Q. 
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Appendix B 

Quantum corrections to conductivity 



B.l Diagrammatic perturbation theory 

Assume the simplest model for the scattering potential a set of randomly distributed im- 
purities with short-range potential, 

U(r) = J2 u ( r ~ r i) ~^J2 6 ( r ~ r ^- 

3=1 3=1 

Here iVj is the number of defects, Uq is the potential amplitude, while V is the system 
volume. The the correlation function of the potential is in the main approximation 

(U(r)U(T>))=n i v*6(r-T>). 

Here rij = iVj/V is the impurity concentration. The following procedure is actually valid 
in the limiting case 

iVj ^oo, V^oo, u ^ , riiul = const . 

Let us assume also zero temperature to avoid inelastic scattering. 

Under these assumptions it is convenient to introduce the resolvent operator 

G{z) = (z-H)- 1 

where z is a complex number, while Ti is the Hamiltonian of the system. Since formally 

e — ri 

where P means the principal value, we can write 

g(e) = ±-T T 5(e-H) = T-^% [Tr G(e ± iO)} . 

113 



114 



APPENDIX B. QUANTUM CORRECTIONS TO CONDUCTIVITY 



The matrix elements of G, 

(r\G(e ± iO)|r'> = £ ^ff^ = G ^ ^ ^ 

a 

are called retarded (+i0) or advanced (— iO) Green functions. They have a simple physical 
meaning of the propagation amplitudes for the states with a given energy, 

Gf = \A(r -> r') "along the time" 

oft 

while 

Gj = -A(r — > r') "against the time" . 
n 

At kp£ ^> 1 a perturbational approach is possible which allows to expand Green functions 
in powers of scattering potential and zeroth-order propagators 

G (0) (z,r-r') = J2 G{0) k ) e ik(r ~ r,) , (B.2) 

k 

where 

G< >(*,k) = (k|#°)(»|k') = -^—5^ . (B.3) 

Z - 6 k 

It is convenient to use k-representation and expand 

^(k)C(k') 



e — e„ ± iO 



The result can be expressed in a diagrammatic form, as shown in Fig. p.l| , panel (a). Each 
solid line corresponds to the non-perturbed propagator ( |B.3|) , each dashed line corresponds 
to the transferred momentum q^, each vertex corresponding to scattering at the defect at 
the point r,- contributes as 

5(k J+l - k, - % ) . 

After n scattering events the incoming wave vector k is changed to k' = k + YTj=i C U- For 
given k and k' one should sum over all keeping the above conservation law. 

The impurity average can be performed by summing over all pairwise contractions [see 
Fig. B.l , panel (b)], each contraction contributes as 

(C/(q i )?7(q J -)) = n^V5 qi _ q . . 

The formal summation of the perturbation series can be performed by introducing the 
irreducible diagrams, E, which cannot be divided into parts by cutting only a sing le G*°>. 
As it can be easily seen, in this way we arrive at the Dyson equation (see Fig. p.2|) , 

(G) = G (0) + G (0)) S (G) . 
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Figure B.l: a - Diagrammatic representation of the perturbations in scattering strength, 
b - Perturbation series for the average Green function. 



+ 




Figure B.2: Illustration of the Dyson equation. 



Its formal solution is 

(G) 



Since this solution sums up an infinite series, we calculate the self-energy E in the lowest 
order in iiiU^, which is just the triangle in the first diagram in Fig. B.l| , panel (b). That 
corresponds to the Born approximation in quantum mechanical scattering. The result is 

E B (z, k) = riiul = Uiul I de . (B.4) 

q Z — £k+q J Z — 6 

Taking z = e + iO the associated retarded self-energy has both real and imaginary part. 
The real part (actually divergent) can be absorbed into renormalization of bare energies 
Ek . The imaginary part 

3 = -irniulgo(e) 

near the Fermi level can be replaced by a constant 

7 = h/2r = hv F /£ 

yielding the average Green function 
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that is equivalent in the coordinate representation to the Green function depending on 
r = r i — Yf as 

This is just damped outgoing spherical wave. 

B.2 Kubo-Greenwood formula 

A time-dependent electric field E(r,t) = — W(r, t) leads to a change in the Hamiltonian 
'Hit) — > T~t — eV(r,t). Let us for simplicity assume that 

V(r,t) =y(r)e _a|t| coswf, s -> +0 . 

Assume that the electric field is concentrated only in the region where the scattering takes 
place and introduce a complete orthonormal set of exact eigenstates ip a (r) indexed by 
continuous variable a. 

To calculate current we need the nonequilibrium density matrix g which is determined 
by von Neumann equation 

The equilibrium state is described by the density matrix 

g(°) = J da f(e a )\a)(a\ , 

where f(e a ) is the Fermi function. The first-order correction is 

mQ W = [H, gV] + [Wi,£ (0) ] . 
Introducing fp a = f{ep) — f(e a ) and V a p = (a\V(r)\fl) we obtain 

ihd t g {1) = -epaffcl - ef Pa V af3 e~ s ^ cos(ut) . 
Its solution for the initial condition g — > at t — > — oo has the form 

Now we can calculate the current, 

j(r,f) =Tr (^j(r)) = J da d(3 j^r) . 

Here 

ipa(r) = ^ [^(r)2ty a (r) - i> a {r)V*r p {v)] , 



B.3. WEAK LOCALIZATION CORRECTIONS 



117 



where 

V = V + (ie/hc) A(r) 
is the covariant derivative. Expressing 

eV aP = — [ d d v j a/3 (r) -E(r) 

£/3a J 

and collecting the contributions to the current as j(r) cos(u>t), we obtain 

j(r) = J d d r' a{r,r') E(r') 
where non-local conductivity a(r, r') is given by the expression 

f/3a 



da dp 



£/3a — hluJ + illS 



+ {uj — >• —a;) 



Here ® means tensor product of two vectors. Having in mind the definitions of the Green 
functions and taking Fourier transform in the quantity r — r' we find that non-local con- 
ductivity at T = can be expressed as 



H s e 2 



3m 2 irV 



} j k i k' i K] c ] l .t(q, u) 



k.k' 



where 



Kw{q,u) = G R (k + ,k' + ,e F + hw)G A (k_, k'_, e F ) . 
Here k±k ± q/2. Note that the function 



(B.7) 

(B.8) 
(B.9) 



k,k' 



has a physical meaning of a density relaxation function. 



B.3 Weak localization corrections 

The function K kk /(q, uj) is a two-particle Green function. We are interested in its impurity 

average, 

/Ckk'(<?,^) = (Kw(q,uj)) . 

Again, we can extract irreducible vertex part, rkk'(q, oS) [see Fig. |B.3| , panel (a)] As a result, 
we rewrite i^-function as 

1CMq>") = G t G t SicM' + G^GiTwi^G^G* . (B.10) 

Here G R ' A denote averaged Green functions which take into account the impurity lines 
which start and finish at the same line. The vertex part F takes into account correlations 
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k+ K k+ k+ k+ 





k. k. k. k _ k >_ 




= Q+ O0+ a + 



Figure B.3: Graphical representation of the disorder average of the product G R G A . The 
thick lines denote averaged Green functions. 

between the quantum mechanical amplitude and the time-reversed one due to scattering 
against the same impurity centers. 

The first term in this equation leads to the Drude formula, while the quantum correc- 
tions are incorporated in the vertex part, V. If the corrections are small, this part can be 
calculated using the main sequence of the diagrams. Two types of diagrams are usually 
discussed - "diffusons" and "cooperons", they leads to the main contributions provided 

M» i- 

The first type is the so-called "ladder" diagrams [see Fig. [B.4| , panel (a)]. Summing the 
"ladder", diagrams we obtain 

rki ' (? ' w) -yRM' (B11) 



where 

C(q, u>) = ^ Yl + MGl (e F ) . (B.12) 

k 

For q, uj — > one can expand £k± ~ £k ± ^(q • v )/2 to obtain 

. f dQ^ 1 9 

C(Q, u = / -:— i : ——/ v ~ 1 + iuJT ~ D( l T 

J Att 1 — tour + zr(q ■ v) 

where D = vp£/3 is the diffusion constant. Thus, at ur,q£ <C 1 

T ™w> u >-^V - iu} + Dq 2- 
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Cooperon 



v 1 , / 



= 0+0 



Figure B.4: On the "diffuson" and "cooperon" contributions. 



One can easily check that the associated density relaxation function has a diffusion pole. 
At the same time, the "diffuson" part does not enter the expression ( |B.7|) at q —>■ 0. for 
conductivity because of asymmetry in k and k'. Actually, this is a consequence of our 
model of point-defect scattering. For more realistic potentials "ladder" diagrams lead to 
renormalization of the elastic mean free time r into the transport relaxation time r tr . 

Now let us turn to the "cooperon" contribution [Fig. [B.4| , panel (b)]. We can obtain 
the result without doing a new work by exploiting time reversal symmetry and assuming [] 

(k'_|G A |k_) = (-k_|G A | -k'_). 

As a result of the invariance of the replacement k <-» — k' 

TkkKcL ^) = r k _ k / +q / 2 ,k'-k+q/2(k + k', u) . 



As a result, at 
one obtains 



Q = k + k' ^ 0, lo^O 

&M = w=dix?- (B13) 

The singularity at k + k' — > is due to backs cattering. The expression (P-13|) is called the 
cooperon because of formal analogy with superconductivity. 

Substituting ( |B.13| ) into general expression for the conductivity ( |B.7|) and replacing 

(hkj) — > kp/2 we get 



hn^ulvp f d 3 k 
3ttt J J2n) 

1 This not true in the presence of magnetic field, or inelastic scattering 



Q 
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Making use of the relation 

d 3 k 



I 



{G R G A f « Anr 3 g(e F ) 



(2tt)* 
we finally get 



M"H-^£— (B.i5) 

To make summation over Q let us use the following trick. Namely, recall that 



' / dte&- DQ2)t 

Jo 



-iu + DO 2 

On the other hand, exp(—DQ 2 t) is the Fourier transform of the Green's function of the 
diffusion equation, 

pdP(r, t)t = PV 2 P(r, t) , P(r, 0) = 5(r) . 
This Green's function has the form 



P(r, t) = — — . exp 



^2 



(AnDty/ 2 1 \ ADt 
Since 

we obtain the formal result, 

2e 2 D /"* max 



irh 



f max rft 

L w- (B16) 



The results depend essentially upon cut-off times, t max , m m- Usually it is chosen £ min = 
t, tmax = Tip- However, there are much more sophisticated ways to treat more realistic sys- 
tems. The usual way is to replace the Green's function P(r, t) of the diffusion equation by 
P(r, £)e _ */ T¥ \ To calculate P(r,t) for a realistic system, a boundary problem is formulated 
for the eigenfunction of the diffusion operator, 

-DV 2 P S = \ S P S , d n P s | vac = 0, P(r, t) | lcad = . 

where subscripts "vac" and "lead" mean interfaces with vacuum and leads, respectively. 
Since 

P(r,t) = ^P s (r)e-^, 



we obtain 



2e 2 P 1 (X ^ T) 1 
Aa = > — 



s ^ 
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Assuming the sample infinite in x, y-direction and L z <C L^ = ^/Dr v we get 



The above mentioned procedure allows a rather simple generalization for the case of mag- 
netic field. It appears sufficient to change the diffusion equation to 

~ D ( V + t fc A ) Ps(r) = AsPs(r) ' 

The boundary condition are also changed to include covariant derivatives, 

d n -> d n + 2i(e/hc)A n , 

as in the Ginzburg-Landau equations for superconductivity. 

Let us apply this procedure to calculate the conductivity of an infinite slab with 
L z <C L v in the perpendicular magnetic filed. To find the eigenvalues we can consider 
a Schrodinger equation for a particle with the mass M = h/2D m/k F £ and charge — 2e 
in the presence of magnetic field. As a result, 



\ s = (s + l/2)(4DeH/h). 



As a result, 



A<7 2d (if) = 



2e 2 DeH 



E 

s=0 



1 AeDH 

r 9 ch 



S + 



where s max w eft/ (2eHl 2 ) 3> 1. In the limit £ <C L v , ^Jch/2eH the result can be expressed 
through the standard Digamma function, 



Acr 2d (tf) 



27T 2 fl 



V 2e#£ 2 / 1 2 2 



eft 



2eHDr ir 



(B.17) 



This expression provides numerical factors. It has two clearly different regimes which cross 
over at 

ft 
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Appendix C 

Derivation of Landauer formula. 



In this derivation I partly follow Ref . [113] . According to quantum mechanics, the expecta- 



tion value of any single-particle operator A can be expressed through the time-dependent 
density matrix, p(t) as 

(A) = Tr {p(t)A} , 
where pit) satisfies the equation of motion 

p = -(i/h)[H, P ]. 

The unperturbed system is described through the Fermi function /(e) as 

Po= da f(e a )\tp a )(ifj a \ . 



Here \ip a ) stand for exact scattering-wave states of the equilibrium system with energy e a . 
The integral is written to emphasize that we deal with continuous spectrum. 

Let us assume that the system is perturbed by a scalar potential with frequency u 
which is turned on adiabatically at t — > — oo, 

cf) oc exp(ut — iut) , v — > . 

We assume that the potential amplitude approaches a constant value in each lead. Then 
we solve the equation for p up to the linear terms in 0. In the limiting case, just as in 
course of derivation of the Kubo formula for extended systems, we get 

(J(r)) = J dr' a(r, r') • E(r') , 

where the local Kubo conductivity tensor is given by the expression 



<7(r,r') = -H f da dp \f(e a )ir6{ep a )+i^V (— ) 

J L e Pa \ e l3aJ 



>j3a 



r)J a ,(r'). (C.l) 
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Here ep a = ep — e a , jp a = f(ep) — f{e a ), f = df/de and V denotes the principal value. 
J/3 Q (r) is the matrix element of the current operator between that exact eigenstates, 

3 Pa(r) = -— [^(r)P^(r) - ^(r)P*^(r)] , 

where 

V = V - (ie/hc)A(r) 

denotes the gauge-invariant derivative. It follows that in the absence of magnetic field the 
principle value vanishes because to time-reversal invariance. 

The current densities should be integrated over the surfaces of the leads. As a result, 
the conductance is given by the expression 

G mn = / dS m I dS n • &(t, t') ■ n 



Here and n are unit vectors normal to the cross-section. 

In the absence of magnetic field one can immediately integrate Eq. (|C.1|) over energy 
and take the limit T —>■ 0. Then the result becomes a sum over discrete states at ep of the 
form 

^X)W"0Mn), (C2) 

where 

Ip a (m) = J dS m -J/3 a (r). (C.3) 

The normalization factor h~ 2 arises from the integration over the energy. 

Since the energies are fixed, the states are characterized by a mode index a, b, and a 
lead index p, q denoting the lead and mode which contains an incoming wave from infinity. 

In lead m, 

4> a (r) = ^(r) + ^t„ p 0+(r), 

cs 
ds 

Here t sc ^ ap is the transmission amplitude for an incident wave in lead p and mode a to 
scatter to lead s and mode c. $j"(r) are the wavefunctions of the infinite perfect leads 
consisting of a longitudinal plane wave traveling away from the sample multiplied by the 
transverse wavefunction for the mode c. If we are interested in Ip a (rn) we keep onpy the 
items which involve the states existing in mth lead, 



c 
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Thus, 



mq \ 



r m ) d x 



mp\ 



mc<—apYmc 



( r m) 



md*—bqrrnd 



r m ) 



Since transverse modes are orthogonal and normalized to a unit flux we get after integra- 
tion over S m , 



In a similar way, 



4/3 H 



&pq&ab&am ^ ^ tinc^bqt 



mc^ap 



fipq^ab^an ^nd^-bq^ 



ndf—ap 



Substituting these expressions into the Kubo formula ( |C.1| ) and using unitarity of sctat- 
tering amplitides and time reversobolity, after rather long algebra we get 



p 2 2 

G =—Trtt ] = —T 
h h 
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Appendix D 

Coulomb blockade in a 
non-stationary case 



Then we can use Kirchoff 's second law for the 2 loops in Fig. |4.5| , 

V e -^ + ^-V 9 = (D.l) 

L'e ^g 

V c -^ + ^-V 9 = 0. (D.2) 

^g 

Combining equations (|4.11|) , (p.l|) and ( p.2|) , we find all the charges, 

Q e = ^ [C c (V e - V c ) + C g V e - C g V g - ne] ; (D.3) 
Q c = ^[-(C e + C g )(V e -V c ) + C g V e -C g V g -ne]; (D.4) 
Q 9 = § [C e (V e - V c ) - (C e + C C )(K-)V 9 - H . (D.5) 



Now let us suppose that 1 electron tunnels from emitter to the grain, then the number of 
excess electrons is changed, 

n — ► n + 1 , Q e — > Q e - e . 

This situation is non equilibrium, so a current must flow through the external circuit. The 
net change of charge with respect to the equilibrium value for n + 1 excess electrons given 
by Eq. ©is 

SQ e = -e^r = -e - 5Q C - 5Q g = -e + e ^ c . 

We observe that the charge transferred through the voltage source V e is ea e e where a ee = 
(C c + C g )/C. Looking at Eq. ( p.4| ) we find that the charge transferred through the voltage 
source V c is a ec e where a ec = C c /C. In a similar way we can consider the jump of 
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the electron from the grain to collector. In this case we have the transferred charges 
ea ce , a ce = C e /C through the source V e and ea cc , a cc = (C e + C g )/C through the source 

Vc ' 

Let us define the tunneling currents through emitter and collector junctions as I e / C (t), 
respectively. Then the current in the emitter branch is 

Ii(t) = a ee I e (t) + a ce I c (t) 

= %^h(t) + ^I c (t). 

In a similar way, the current in the collector lead is 

I r (t) = a ec I e (t) + a cc I c (t) 

= -jjh(t)+ c 9 Ic(t). 

Usually, C g <C C e , C c , and I t fa I r m I(t), where 

C C 

I(t) = a e I e (t) + a c I c (t) , a e = , a c = . (D.6) 

The partial currents can be expressed through the probability p n (t) to find n excess elec- 
trons at the grain and the tunneling rates T^ u (n), 

U*) = e - r^ M ] . (D.7) 

n 

The probability p n (t) should be calculated from the Master equation, 

^ = Pn-lK-l + Pn + lK+l ~ (^^ + Pn (t) . (D.8) 

Here 

= r e ^ g (n - 1) + r^ g (n - 1) ; (d.9) 

T n n+1 = r^ e (n + l) + r^ c (n+l). (D.10) 
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